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The fundamental properties of the lepton sector include the neutrino masses and flavor mixings. Both are 
difficult to observe because of the extremely small neutrino masses and neutrino-matter cross sections. In these 
lectures, we focus on the basic concepts for the determination of neutrino properties. We introduce neutrino 
oscillations as standard mechanism for neutrino flavor changes, and we discuss methods to measure the neutrino 
mass. Furthermore, we illustrate how precision measurements in neutrino oscillations will be performed in the 
future, and may even open a window to new physics properties, such as motivated by LHC physics. Finally, we 
discuss some applications of neutrinos in astrophysics, such as neutrino oscillations in the Sun. We also illustrate 
how neutrinos from extragalactic cosmic accelerators may be used for the determination of neutrino properties. 



1. Introduction 

Neutrinos are the most abundant known mat- 
ter particles in the universe, in number, exceed- 
ing the constituents of ordinary matter (electrons, 
protons, neutrons) by a factor of ten billion. The 
neutrino fluxes are extremely high. For exam- 
ple, there are about 7 • lO^^cm^^s^^ streaming 
through the Earth from the Sun. Neutrinos are 
naturally produced in the Big Bang, the Earth's 
crust (by uranium and thorium decays), the Sun, 
supernovae, and the Earth's atmosphere by cos- 
mic ray interactions. Furthermore, at the high- 
est energy scales, they are most probably pro- 
duced in cosmic accelerators together with the 
observed cosmic rays. Man-made neutrinos are 
emitted from nuclear fission reactors, as well as it 
is possible to produce artificial neutrino beams. 
Although they are very abundant, it is difficult to 
catch them because of their extremely small cross 
sections. As a consequence, they even escape 
from very dense environments, such as super- 
novae or nuclear power plants. Roughly speaking, 
the number of detected neutrinos scales as 

A^obs a (/) X cr X t X TTlDet ■ (1) 
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The flux (/) is extremely large, the cross section 
cr extremely small. The observation time should 
not be longer than a few years. In order to accu- 
mulate sufficient statistics, one therefore can esti- 
mate the necessary detector sizes from Eq. ([T|) to 
be of order of kilotons. In the Standard Model 
(SM) of elementary particle physics, neutrinos 
are massless particles. However, after a long his- 
tory of disputed alternatives, neutrino oscillations 
have recently been established as leading neutrino 
flavor change mechanism, which implies that at 
least two of the neutrino (mass) states must be 
massive. Massive neutrinos are now one of the 
very few specific hints for physics beyond the SM. 

Although the SM can be easily extended by 
right-handed neutrinos to introduce Dirac mass 
terms, the lightness of the neutrinos then points 
to un-plausibly small Yukawa couplings. A typ- 
ical way out is the introduction of a Majorana 
neutrino mass, which implies that the neutrino 
is its own antiparticle. This alternative has a 
number of interesting implications. First of all, 
this hypothesis can be, in principle, tested in 
neutrinoless double beta (Ozy/3/3) decay, unless it 
is suppressed by intricate parameter constella- 
tions. Second, a Majorana mass term violates 
lepton number, which is an accidental symme- 
try of the SM. Such a lepton (or baryon) num- 
ber violation is, for instance, needed in dynami- 
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cal mechanisms to describe the observed matter- 
antimatter- asymmetry in the universe (baryoge- 
nesis). Maybe neutrinos are somehow involved in 
such a mechanism (leptogenesis (Tj). And third, 
neutrino mass can be interpreted as the low- 
est order perturbation of Beyond the SM (BSM) 
physics in the sense of the effective operator 
picture. In this picture, the famous Weinberg- 
Operator [2] 

■^elr^ oc J (LicT^H) (Lia^H) +h.c. (2) 

(L: lepton doublet, H: Higgs doublet) is the only 
d = 5 operator, and it is the only operator sup- 
pressed by only one power of A, where A is the 
BSM physics scale. It leads to Majorana neu- 
trino masses after electroweak symmetry break- 
ing (EWSB). Using order one couplings and neu- 
trino masses < eV, one can easily show that A in 
Eq. ([5]) points towards the Grand Unified Theory 
(GUT) scale, the scale where the gauge interac- 
tions of the SM are expected to unify. In this 
case, it may be obvious to suspect a connection 
with the mentioned baryogenesis concept in the 
early universe. The fundamental theories leading 
to the operator in Eq. ^ and their implications 
are discussed in the lecture by Z.-z. Xing in this 
series [3] (see also [4] for some aspects). 

In order to describe possible theories for neu- 
trino mass and the connected BSM physics, it is 
mandatory to pin down the properties of the neu- 
trinos. As we shall see, most of the observables, 
such as the mass (squared) splittings and mixing 
parameters, are accessible in neutrino oscillations. 
In particular, the observation of leptonic CP vio- 
lation may support the argument that neutrinos 
are involved in leptogenesis. Neutrino oscillations 
can be used to fix essentially all expected (guaran- 
teed) observables except from the absolute neu- 
trino mass scale and possible Majorana phases. 
The direct measurement of the neutrino mass typ- 
ically involves the determination of the endpoint 
in the electron spectrum coming from tritium de- 
cays. The direct test of the nature of neutrino 
mass involves the test of Oz//3/3 decay. In addition, 
upper bounds on the neutrino mass scale can be 
obtained from cosmology. Apart from neutrino 
flavor mixing and neutrino masses, other proper- 



ties of the neutrinos may point towards the na- 
ture of neutrino mass, such as the electromagnetic 
dipole moment, or to new physics effects, such as 
neutrino lifetime. Therefore, it is worthwhile to 
test and constrain these as well. Especially neu- 
trino from extragalactic cosmic accelerators may 
allow for probes in different baseline and energy 
regimes, where new physics effects may be present 
which are otherwise not observable. 

In these lectures titled "neutrino phenomenol- 
ogy" , we focus on the determination of the neu- 
trino properties. In Sec. [21 we discuss the current 
understanding of neutrino oscillations. Then in 
Sec. m we show future perspectives for neutrino 
oscillation experiments, with special emphasis on 
the discovery of leptonic CP violation. In this 
section, we add matter effects in (constant) Earth 
matter to the framework of neutrino oscillations. 
As the next step, we illustrate how neutrino oscil- 
lations work in varying matter density in the Sun 
in Sec. |4l and show how these can be used to test 
the solar neutrino mixing angle. In Sec. [5] we dis- 
cuss the possibility to use neutrinos from cosmic 
accelerators for the test of neutrino properties, 
such as neutrino lifetime. Finally, we summarize 
the approaches to test the absolute neutrino mass 
scale in Sec. [H 

2. Neutrino oscillation framework 

Here we present the current understanding of 
neutrino oscillations. We start with a short his- 
torical perspective, then introduce neutrino oscil- 
lations in vacuum in the standard quantum me- 
chanical treatment, and comment on leptonic CP 
violation. Then we derive the two-flavor limit, 
and demonstrate that the general three-flavor 
case can be reduced to two-flavor sub-sectors for 
our current understanding. 

2.1. Historical perspective 

Historically, neutrinos from the Sun were ob- 
served for many decades by Raymond Davis, Jr. 
and collaborators since the 1970s in the Home- 
stake experiment 5 . On the other hand, solar 
models based on the nuclear fusion chains in the 
Sun, such as by John N. Bahcall and collabora- 
tors (see, e.g., Ref. [5] for a recent discussion). 
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predicted much higher (electron neutrino) fluxes 
based on the normalization to the solar luminos- 
ity. This solar neutrino anomaly (see Ref. [7] 
for an early reference) can be plausibly described 
by neutrino flavor changes. The most important 
results in solar neutrino physics in this context 
is probably the SNO (Sudbury Neutrino Obser- 
vatory) neutral current measurement [8^, which 
confirmed the flux predictions of the solar model, 
and therefore the flavor changes of the neutrinos. 
Since the neutral current interactions are equally 
sensitive to all (active) neutrino flavors, the ra- 
tio between charged current interactions of elec- 
tron neutrinos and neutral current interactions 
directly determines the fraction of neutrinos still 
found in the original state. Apart from neutri- 
nos from the Sun, neutrinos are abundantly pro- 
duced in the Earth's atmosphere. Cosmic rays, 
which could be protons or heavier nuclei, hit the 
Earth's atmosphere to produce showers including 
charged pions. These charged pions decay with 
the decay chain 

7r+ ^+ + i^^ , 

/i+ ^ e+ -I- -f i^e , (3) 

as illustrated for 7r+ here, leading to a flux of elec- 
tron and muon neutrinos and antineutrinos. The 
most compelling evidence for neutrino oscillations 
as leading flavor change mechanism is probably 
coming from the observation of atmospheric neu- 
trinos in the Super-Kamiokande experiment [5]. 
This experiment has used the directional depen- 
dence of the incoming neutrinos, which deter- 
mines the path length traveled through the Earth, 
to infer the neutrino oscillation parameters. In 
fact, it has turned out that solar and atmospheric 
neutrino flavor changes can be described by dif- 
ferent sets of oscillation parameters. The cor- 
responding oscillations can be described as sub- 
sectors of the general three-flavor case. 

In this section, we choose a deductive rather 
than historical presentation. We show that the 
existing measurements can be described by lead- 
ing two-flavor sub-sectors, which our understand- 
ing is based upon. Although such a presenta- 
tion suggests that neutrino oscillations have been 
taken for granted as flavor change mechanism, 
they have been established by considering a num- 



ber of alternatives, such as neutrino decay and 
decoherence, and a number of anomalies. For a 
more refined presentation of the current picture, 
see Ref. [TU], and for a more detailed presenta- 
tion of the neutrino oscillation phenomenology, 
see, e.g., Ref. [11 and references therein. In the 
next section, we will then look beyond these two- 
flavor sectors and introduce three-flavor effects. 
Although some of the recent experiments are al- 
ready sensitive to three-flavor effects, we perform 
this splitting for conceptual reasons. 

2.2. Neutrino oscillations in vacuum 

In a quantum mechanical picture, the eigen- 
states of the weak charged current interactions 
Iva) do not correspond to the mass eigenstates 
|z^fc), i.e., the eigenstates of the free Hamiltonian 
H with the neutrino energy eigenvalues Ek = 
y7? + mf . The flavor and mass eigenstates are 
connected by a unitary {N + S) x (N + S) matrix 

N+S 
k=l 

where N is the number of active and S is the num- 
ber of sterile (i.e., not weakly interacting) neu- 
trino mass eigenstatesH Greek indices denote fla- 
vor eigenstates, Latin indices denote mass eigen- 
states. The unitary mixing matrix U is, for iV = 3 
and S = 0, also called C/pmnSj or Pontecorvo- 
Maki-Nakagawa-Sakata matrix. The flavor and 
mass eigenstates, respectively, are both assumed 
to form a basis. Applying the time evolution of 
the mass eigenstates in vacuum 

\vk{t)) = ex.Y>{-iEkt) \vk) (5) 

to Eq. Q and using the unitarity of the mixing 
matrix, we obtain the vacuum transition ampli- 
tude 

Au^^up = = {vpWait)) (6) 

N+S 

= ^-kUpk exp{-tEkt) . 

fc=i 

^The contribution of sterile states is strongly constrained. 
However, we keep them in the derivation to demonstrate 
where the calculation works independently of the number 
of participating flavors. 
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Note that this equation can be equivalently ob- 
tained from the Schrodinger equation in matrix 
form 



ity in vacuum 



at 



El 
E2 



(7) 



with "Hf the Hamiltonian in flavor space and ^ 
the flavor state vector. This form is often used 
if the Hamiltonian is explicitely time-dependent, 
such as it might be for the matter effects discussed 
later. For the transition probabihty, we have from 
Eq. O 



Pc 



KpAo^p = (8) 



N+S 



k,j=l 



The quantity J^^ is also known as quartic re- 
phasing invariant [12\ , which characterizes the in- 
formation in the mixing matrix independent of a 
possible phase re-definition of the charged lepton 
and neutrino fields. The standard derivation of 
the oscillation formula relies on the approxima- 
tions for ultra-relativistic neutrinos 



El, 



E 



2E 



(9) 



which imply that the different mass eigenstates 
have different energies. This seems to be con- 
tradictory to a neutrino produced with a certain 
energy as a superposition of mass eigenstates. 
Therefore, in any more refined calculation, as- 
sumptions on the energy and momentum widths 
have to be made. The simplest such method in- 
cludes the assumption of wave packets, see, e.g., 
Ref. [13], which leads to the same oscillation for- 
mula we will obtain - as long as there is enough 
wave packet overlap among the different mass 
eigenstates. From Eq. ([5|) using Eq. © and the 
definition Am^^- = m^, — m|, we find after some 
transformations the neutrino oscillation probabil- 



P, 



= - 4 ^ Re J^/ sin2 ( 

k>3 \ 



AE 



k>j 



CP conserving 

-^sin(^]. (10) 



2E 



CP violating 

The quantity L, the distance between source and 
detector, is often called baseline, and the func- 
tional dependence J-{L, E) on L and E is of- 
ten called spectral dependence; in vacuum, it is 
just F{L,E) — L/E, as one can read off from 
Eq. (flOl) . The evidence for neutrino oscillations 
as leading flavor change effect comes from this 
particular spectral dependence, i.e., the flavor 
change effect as a function of energy and baseline. 
There is no sensitivity of neutrino oscillations to 
the absolute mass scale, but the mass (squared) 
splittings and the ordering of the masses is, in 
principle, determined by Eq. ([TU]). For N -\- S 
flavors, there are N + S ~ I independent mass 
squared splittings. 

2.3. On leptonic CP violation 

A very important quantity of interest is the CP 
symmetry, where "CP" stands for charge-parity. 
It is basically a symmetry between the behavior 
of particles and anti-particles taking into account 
the peculiarities of the electroweak framework (in 
particular, the V — A interactions only coupling 
to left-handed particles and right-handed anti- 
particles). In the context of our CP asymmetric 
universe in which we do not find sufficient an- 
timatter to justify the CP symmetry, the ques- 
tion of the source of the violation of this sym- 
metry, in short, CP violation, is probably one of 
the most interesting ones in particle physics. The 
CP conserving part in Eq. (fTO|) is the same for 
neutrinos and anti-neutrinos. The CP violating 
part changes sign for anti-neutrinos (for which 
the mixing matrix effectively has to be complex 
conjugated). It is only present if the mixing ma- 
trix has complex phases (apart from possible Ma- 
jorana phases) invariant under a re-definition of 
the lepton fields. Note that the CP violating 
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part oscillates with the double frequency com- 
pared to the CP conserving part. In addition, 
note that ImJ^^ is, up to a sign depending on 
the indices, equivalent to the so-called Jarlskog 
invariant J [141 , which is frequently used for the 
quantification of CP violation. Another quantity, 
which has been used historically, is the CP asym- 
metry 



P. 



a-cp 



af3 



Pap ' 



(11) 



where refers to Pp^^p^. For modern statis- 
tical simulations, this quantity is not very repre- 
sentative, mainly because of matter effects, which 
violate the CP symmetry extrinsically by the ab- 
sence of antimatter in the Earth, as we will dis- 
cuss in the next section. However, one easily 



finds from Eq. (|T0]) and Jj 



U, 



ak\^ \Uaj\^ that 



acp = for a — (3, i.e., one needs to observe the 
transition among flavors to access CP violation. 
If the oscillations are averaged out, such as by a 
poor energy resolution of the detector, we have 
for the oscillatory terms in Eq. ([T0|) 



2 / ^^hL 

AE 



L/E 



sm 



2E 



_ 1 
" 2 ' 

= 0. 



(12) 



L/E 



This implies that the measurement of CP viola- 
tion requires in addition the observation of the 
spectral dependence. 

2.4. TviTO-flavor limit 

In order to illustrate neutrino oscillations, it is 
useful to consider the two-flavor limit, i.e., N = 2 
and 5 = 0. From the simple two-flavor mixing 
matrix 



U = 



cos 9 sin 6 
— sin 9 cos 9 



(13) 



parameterized by only one mixing angle 9, we di- 
rectly obtain from Eq. pop the transition proba- 
bilities for the two flavor states \va) and sep- 
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Figure 1. Disappearance probability Paa as a 
function of x = L/E in the two-flavor limit. The 
different curves represent different values of the 
mixing angle for illustration. The shaded re- 
gion corresponds to the appearance probability 

Pa/3 = 1 - Paa (for 9 = 14°). 



arated by only one mass squared splitting Ato^: 



Pa 



1 - sin^ 29 SIT? 

BIT? 29 BIT? 



\ 4E 
An?V 



AE 



(14) 



The first probability is often called disappearance 
probability or survival probability (because the fla- 
vor Va disappears or survives, depending on the 
point of view) , and the second probability is often 
called appearance probability (because the flavor 
lyp appears). These probabilities are visualized 
in Fig. [T] The mixing angle sin^ 29 can be inter- 
preted as the oscillation amplitude, whereas the 
mass squared splitting Am^ can be interpreted 
as oscillation frequency, and its inverse is propor- 
tional to the oscillation length A = (47rP)/ATO^. 
The two probabilities add up to one, which is 
a consequence of the unitarity of U. Note that 
even if one introduces an additional CP phase in 
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Figure 2. Neutrino mass eigenstates for normal 
and inverted mass ordering (not to scale). 



Eq. (|13|) . the quartic invariant J^^ in Eq. (fTO|) 
cannot become complex for two flavors, which 
means that there will be no CP violation observ- 
able in two-flavor oscillations. For the same rea- 
son, i.e., the presence of CP violation in flavor 
mixing, Kobayashi and Maskawa postulated three 
flavors in the quark sector, for which they received 
the Nobel prize 2008. 

2.5. Three-flavor case 

The current standard assumptions for neutrino 
oscillations include three active and no sterile 
neutrino flavors, i.e., N — 3 and S ~ 0. In 
this case, the possible mass spectra are illus- 
trated in Fig. [2] for a normal {Ani^^ > 0) and 
inverted (Am^i < 0) mass ordering. The split- 
tings between the mass eigenstates are deter- 
mined by Am^x = rn^ — ml, also called so- 
lar mass splitting, and Am§]^ = m§ — to^, also 
called atmospheric mass splitting (Arri^2 is given 
by Am|j^ — Am2i). Since the upper bound for the 
neutrino masses is known to be of the order eV, 
as we will discuss later, the mass spectrum can 
be close to this bound, called degenerate spec- 
trum, or close to zero, called hierarchical spec- 
trum. Since neutrino oscillations are not sensitive 
to this feature, the terms "normal/inverted or- 
dering" and "normal/inverted hierarchy" are of- 
ten used equivalently. The mass ordering and the 
type of the spectrum are characteristic for neu- 
trino mass models, see Ref. \I5\. For example, 
the structure of the neutrino mass matrix is qual- 
itatively different for the normal and inverted or- 
dering in case of a hierarchical spectrum. 



The three-flavor mixing matrix is, apart from 
possible Majorana phases not relevant for neu- 
trino oscillations, typically parameterized as |16] 



C^PMNS — 




Atmospheric mixing 

ci3 si3e-'*cF 





-si3 e 



iScF 





Cl3 



Reactor mixing 

Cl2 Sl2 
-Si2 Cl2 
1 



(15) 



Solar mixing 



where dj — cos 6ij and Sij = sin Oij . This implies 
that the neutrino mixing can be parameterized by 
three mixing angles 623, ^13, and 9i2, which are, 
for historic reasons, often called atmospheric mix- 
ing angle, reactor mixing angle, and solar mixing 
angle, respectively. In addition, there is one phase 
Sqp, which leads for S^p ^ {0, tt} to CP violation, 
c/., Eq. PU)) . Note that in this parameterization, 
exp(zi5cp) is multiplied by sin0i3, which means 
that a non-zero value of ^13 is required for any 
measurement of Sep. 

Together with the two independent mass 
squared differences Am^i and Ani^^, we have six 
neutrino oscillation parameters. The mixing ma- 
trix can be fully described by mixing angles in the 
parameter ranges 6ij e [0,7r/2] and 6cp G [0, 27r[ 
(see, e.g., Ref. [H]). If the neutrinos are Majo- 
rana particles, the mixing matrix should be re- 
placed by 



PMNS 



PMNS 



X diag(l,e*",e*''), (16) 



because these additional phases cannot be ab- 
sorbed in a re-definition of the Majorana fields. 
These phases with the physical parameter ranges 
[0, tt] enter the description of OvfUf] decay, but not 
into neutrino oscillations. They are called Majo- 
rana phases, whereas i5cp is often referred to as 
Dirac CP phase (meaning that it is also present 
for Dirac neutrinos). Note that the parameteri- 
zation in Eq. (jlSp is somehow arbitrary. It only 



7 



Parameter 



Best-fit Degrees 



2cr 



3(7 



Am|i [lO-^eV^] 


' •00-0.20 




7.25-8.11 


7.05-8.34 


lAmiil [10-3eV2] 


^■^^-0.11 




2.18-2.64 


2.07-2.75 


sin^ 6*12 


u.ou^_o.oi6 


33° 


0.27-0.35 


0.25-0.37 


sin^ 6*23 




45° 


0.39-0.63 


0.36-0.67 


sin^ 6*13 


^•^J^-o.oii 


6° 


< 0.040 


< 0.056 


^CP 


Currently no 


information 





Table 1 

Current best-fit values with la errors, best-fit values in degrees (angles only), and 2(t and 3(t intervals 
(1 d.o.f.) for the three-flavor neutrino oscillation parameters from global data; adopted from Ref. fT7|. 



makes sense in combination or comparison with 
the quark sector, where Vckm describes the rota- 
tion between the up- and down-type quark states 
using the same parameterization. For example, 
one may test a possible connection between the 
quark and lepton sectors, such as by a unifying 
theory, or obtain hints for the generation of the 
flavor structure, which may or may not have the 
same origin in both sectors. 

The current knowledge on the three-flavor neu- 
trino oscillation parameters is summarized in Ta- 
ble [TJ We can read off two qualitative observa- 
tions from this table, which will be relevant for 
our analytical treatment: 

1. One of these mass squared differences is 



much smaller than the other two: Ar 



Am^i 



Ar 



This leads to a hierarchy of 



the neutrino mass splittings, as illustrated 
in Fig. El 

2. Two of the mixing angles, 623 and 612, are 
very large, whereas one mixing angle ^13 is 
small - at most of the size of the Cabibbo 
angle 9c in the quark sectorial 

There might be even maximal mixing 623 = 7r/4, 
which could point (possibly with a vanishing O13) 
towards a fundamental symmetry between and 



^Recently, a 1.6a claim for 613 > has been made from 
the global analysis of all oscillation data 19j. However, 
this claim depends on details of the analysis and may very 
well come from statistical fluctuations, see Refs. I20I21| for 
a more detailed discussion. 



i^T0 For 012, maximal mixing is excluded at 
more than 5a. However, the mixing angles are 
compatible with the so-called tri-bimaximal mix- 
ing [52], where sin^ 612 = 1/3, sin^ ^23 = 1/2, and 
sin^ 6*13 = leading to a specific form of the neu- 
trino mass matrix, which has been motivated by 
a large class of models. From Table [l] one can 
also read off the primary remaining quantities of 
interest for future experiments: 

• The value of ^13, and if it is non-zero. 



713 

• The sign of Am§2, 
neutrino masses. 



the ordering of the 



• The value of Sqp (only if ^13 > 0), and if it 
is violating the CP symmetry, i.e., if dcp ^ 
{0,^}. 

• The exact value of ^23, in particular, if max- 
imal mixing 623 — -k /A can be excluded, and 
if O23 > 7r/4 or < 7r/4, the 823 octant. 

Above we have mentioned that there is currently 
no evidence for additional sterile neutrino species 
or other new physics contributing to neutrino os- 
cillations in a leading role. For example, the evi- 
dence for active-sterile neutrino oscillations from 
the LSND experiment [53] has been ruled out [53] . 
This, however, does not mean that it is not inter- 
esting to look for sub-leading new physics effects 
in neutrino oscillations, since particular classes of 
effects might be primarily visible there. 



■^Maximal mixing is, from the oscillation point of view, 
illustrated in Fig. [l] In this case, the two-flavor survival 
probability may even vanish at certain L and E. 
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2.6. Two-flavor sub-sectors 

In the following, let us use the qualitative 
knowledge on the neutrino oscillation parameters 
in order to reconstruct the different two-flavor 
sub-sectors which have lead to the current knowl- 
edge. This is not meant to be a complete review, 
but only a short discussion to give the reader 
some idea of the relevant measurements. For the 
sake of simplicity, let us first of all assume that 
Am2i ^ ^nd f/pMNS is real. Then we have 

from Eq. 1^0^ 

Pa/3 = <5„/3 - 4 ( Jgf + Jgf ) Sin^ A31 

- 4J2fsin2A2i, (17) 

where = Am|^ L / (4i?) . We can now choose 
one of the following two oscillation frequencies: 

The atmospheric oscillation frequency or 

A31 ~ 7r/2. This necessarily leads to 
A21 ^ 1, i.e., the second oscillatory part 
in Eq. (1171) is very small. 

The solar oscillation frequency or A21 ~ 

7r/2. This necessarily leads to A31 ^ 1, 

meaning that the first oscillatory part in 
Eq. p7)) averages out; see also Eq. (fT2|) . 

Note that we "choose" an oscillation frequency by 
the neutrino energy E, determined by the neu- 
trino source, and the baseline L, determined by 
the experimental configuration. In addition, note 
that the names "solar" and "atmospheric" fre- 
quency or mass squared splitting (if referring to 
the corresponding Am^) has again historical rea- 
sons, as we shall see below. In the limit ^13 — )■ 0, 
the mixing matrix in Eq. (1151) simplifies to 

(C12 S12 \ 

-S12C23 C12C23 S23 . (18) 
S12 S23 -C12 S23 C23 / 

Using Eqs. (fT7)) and ((T5)) . one easily obtains the 
leading order two-flavor oscillation probabilities 
of the following experiment classes, which have 
been carried out so far: 

Atmospheric experiments use the neutrinos 
produced (mainly) by pions as secondaries of cos- 
mic ray interactions in the Earth's atmosphere. 



Charged pion decays produce fluxes of elec- 
tron and muon neutrinos (and anti- neutrinos), 
see Eq. (|31). The detectors, such as Super- 
Kamiokande [9 , can detect electron or muon neu- 
trinos, which means that the following four oscil- 
lation probabilities are interesting: 

^ ee — ^7 efi — ^ — ^ 

P^^ ~ 1 - sin^ 2023 sin^ A31 . (19) 

Obviously, atmospheric neutrino oscillations can, 
to leading order, be described by the two-fiavor 
limit with the parameters 623 and Am|]^ (in fact, 
the neutrinos change flavor into Vr, which are in- 
visible to the detector). Therefore, these oscilla- 
tion parameters are often called atmospheric pa- 
rameters. 

Solar experiments historically detect the 
neutrinos produced by fusion reactions in the 
Sun, which cannot be described by the vacuum 
oscillation framework we have introduced so far. 
However, we can describe a very long baseline 
reactor experiment using multiple nuclear power 
plants in Japan as neutrino sources: the Kam- 
LAND experiment ( "Kamioka Liquid scintillator 
Anti-Neutrino Detector") [5S]. Since nuclear re- 
actors produce only (by beta decays), which 
might be detected by inverse beta decays, the 
applicable oscillation probability from Eqs. (fTT]) 
and dH]) is 

Pee ~ 1 - sin^ 26I12 sin^ A21 . (20) 

The parameters measured in this experiment are 
the solar parameters, and the probability again 
corresponds to the two-flavor limit. In fact, here 
the oscillate into a superposition of about equal 
amounts of Dfj_ and Dr (with their ratio determined 
by 023) ■ 

Reactor experiments for O^s- Relaxing the 
condition 613 — > 0, i.e., using Eq. (fT5|) instead of 
Eq. (dH), one obtains different from Eq. 

Pee ~ 1 - sin^ 26*13 sin^ A31 

— cos 6I13 sin^ 26I12 sin^ A21 . (21) 

Choosing a much shorter baseline than for the 
solar reactor experiments above, one selects the 
atmospheric oscillation frequency in order to ob- 
tain Pgg ~ 1 — sin^ 2013 sin^ A31, with small devi- 
ations to be interpreted as a signal for a non-zero 
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Figure 3. Reactor electron antineutrino disap- 
pearance probability, illustrated for different sets 
of parameters and perfect resolution. 



^13. Therefore, ^13 is also often called reactor an- 
gle. An example for a corresponding experiment 
has been the CHOOZ experiment (named after its 
site) . The interplay between atmospheric and 
solar oscillation frequency is illustrated in Fig. [31 
The thick solid curve corresponds to 613 = 0, i.e., 
Eq. (PO)) . In this case, the oscillation dip is at 
about 60 km in the gray-shaded solar oscillation 
window. Small changes in Arri^i will move this 
dip, as illustrated by the dashed curve by increas- 
ing Am2i by 20%, which means that such an ex- 
periment is very sensitive to Arn^i. If O13 > 
(thin solid curve) the faster atmospheric oscilla- 
tion will be superimposed, c/., Eq. ([2T|) . leading 
to sensitivity to ^13 in the atmospheric oscillation 
window at about 1-2 km. At the longer baselines, 
the atmospheric oscillations can in practice not be 
resolved and are averaged out. 

Conventional neutrino beams. In this case, 
the neutrinos are produced (mainly) by pion de- 
cays such as in the atmosphere, but using a man- 
made neutrino source. They are detected as 
electron and muon flavors, such as in the cur- 



rently running MINOS experiment ("Main Injec- 
tor Neutrino Oscillation Search", Fermilab) [27] . 
or as tau neutrinos, such as in the OPERA ex- 
periment ("Oscillation Project with Emulsion- 
tRacking Apparatus") in the CNGS ("CERN to 
Gran Sasso" ) beam [28] . The probabilities of in- 
terest are the same as in Eq. ([T9|) (except for 
OPERA, where P^r — 1 — -P^^), for example, 
MINOS has provided an improved measurement 
of ArTigj^. However, with the increasing statistics 
of such experiments, corrections from P^e — 0, as 
in Eq. (|19|) . can be measured. As we will demon- 
strate later, these corrections are not only a mea- 
sure of 9i3 , but also contain the information nec- 
essary to extract CP violation. 

In summary, neutrino oscillations have so far 
mostly been derived from two-flavor sub-sectors 
of the general three-flavor framework. Depend- 
ing on the experiment, the two- flavor probabil- 
ities in Eq. ([H)) are described by different sets 
of parameters, such as {Am|]^,023} (atmospheric 
parameters) for the atmospheric experiments, 
{ATO21, ^12} (solar parameters) for the long base- 
line reactor experiments, and {Am§i, 6*13} for the 
short baseline reactor experiments; see also Fig.[T] 
for typical values of the mixing angles (c/., Ta- 
ble [TJ. Especially the measurement of Sqp will 
be a direct test of the three-flavorness of neutrino 
oscillations, as we shall see in the next section. 

3. Future precision oscillation physics 

In this section, we discuss neutrino oscillations 
beyond the two-flavor sub-sector measurements, 
which have lead to the current knowledge. We 
introduce matter effects in Earth matter to neu- 
trino oscillations, and we show how three-flavor 
effects can be accessed. Furthermore, we intro- 
duce future experiment classes and discuss their 
simulation. Finally, in the era of precision neu- 
trino physics, we also show examples of inter- 
esting new physics effects, and how they can be 
tested. 

3.1. Matter effects in neutrino oscillations 

In order to discuss future precision neutrino os- 
cillation physics, we need another key ingredient 
of neutrino oscillations, which is the matter ef- 
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feet [29130131] , This effect implies that coherent 
forward scattering in matter by charged current 
and neutral current interactions affects neutrino 
oscillations. Neutral current interactions occur 
for all (active) flavors, leading to an overall phase 
which can be subtracted, whereas charged current 
interactions are only possible for electron neutri- 
nos (or anti- neutrinos). The reason for this asym- 
metry is that ordinary matter consists of elec- 
trons, protons, and neutrons, whereas there are 
no muons and tauons required as SU(2) counter- 
parts of the Vfj^ and Vr for charged current inter- 
actions. This leads to an effective net potential 
Aqc on the electron flavor, which can in flavor 
space be described as a term adding to Eq. ([7]) in 
the electron flavoiQ: 

-Hp = U \ I ijt 

I 





2E 



Vcc 

ODD 
ODD 



(22) 



Here Vcc = ±V2Gf rig is the matter poten- 
tial with rig ~ Yep/miy ~ p/{2mN) the elec- 
tron density in matter. The quantity de- 
scribes the number of electrons per nucleon with 
the nucleon mass m^r. Furthermore, the sign 
in Vcc is positive (negative) for neutrinos (anti- 
neutrinos). The evaluation of Eq. ([7^) with 
Eq. (1^^ is straightforward if the Hamiltonian is 
not explicitely time-dependent. In this case, one 
simply re-diagonalizes Eq. ((22)) in order to obtain 
the mixing matrix and mass eigenstates in mat- 
ter (see, e.g., Ref. [32]) ■ This approach, is, in 
fact, often used in numerical calculations, where 
one typically evaluates the Hamiltonian in layers 
of constant matter density. Analytical computa- 
tions, on the other hand, are relatively simple in 
the two-flavor limit in both constant and (slowly 
enough) varying matter densities. In the Sun and 
in supernovae, where matter effects are especially 
important because of the extremely high electron 
densities, one has to deal with varying matter 
densities. Here we first concentrate on the simpler 
case of Earth matter. As long as the baseline does 



^Note that, compared to Eq. 10, we have used already 
the ultra-relativistic approximation Eq. l(9]l here, and we 
have subtracted an overall phase factor. 



not cross the Earth's core, i.e., L < 10 700km, us- 
ing a constant matter density is a good first order 
approximation. In the two- flavor limit, we obtain 
from Eq. (j22p by multiplying out and subtracting 
a global phasqfl 

_ J_ / - Am^ cos 26 + Acc Am^ sin 261 
4E \ Am^ sin 261 Am^ cos 261 - Acc 

(23) 



with 

Acc = 2EVcc^ ±2\/2 EGpUe. 
Compared to vacuum, where 

1 / -Arn2cos26' /\m^ s\n2e 



(24) 



^ Ae[^ ATO2sin26' Am^ cos 261 ) ' "^^^^ 

we can use the same form in matter using effective 
parameters Ajfi^ and 6 

nj _ _}_( -Arh'^ cos2e Arh'^ sm2e \ , , 
~ 4^; V sin 29 Am^ cos 26i j ' ^ ' 

leading to the same form of the oscillation prob- 
abilities Eq. (dH): 



Paa = 1 - sin^ 20 sin^ 



Am^L 



Pap 



sin^ 20 sin^ 



AE 



(27) 



From the comparison between Eq. (|26p and 
Eq. (|23|) . one can easily show that the parame- 
ter mapping is 



Am 



^^2 



Am^ ^ , sin 261 



sin 26* 



cos26l-i) + sir? 20, (28) 



with 



jY - ^^2^/2EGFne 



Aw? 



Aw? 



(29) 



These formulas are useful to demonstrate a 
number of important consequences. First of all. 



^Adding or subtracting to "Hp a. quantity proportional to 
the unit matrix leaves the oscillation physics unchanged. 
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in the limit Am^ L / (AE) ^ 1, the oscihating term 
in Eq. (|27p can be expanded, and one immediately 
can read off these equations that the ^ parameters 
cancel, and the vacuum probabilities are recov- 
ered. This means that long enough baselines are 
relevant to observe matter effects. Second, one 
can read off from Eq. (|28p that the oscillation an- 
gle (or amplitude) becomes resonantly enhanced 
for cos 26* — A. The corresponding resonance en- 
ergy is given by 

Ercs [GeVl - 13 200 cos 26 , / J . (30) 

P [g/cni3] 

For p — 3.4g/cm^ (the average density for L — 
4000 km), Am^ = Aml^ from Table [2 and 
— 6*13 ~ 0, this evaluates to a resonance energy 
of E'res — 9.3 GeV. Therefore, relatively high neu- 
trino energies are required for substantial Earth 
matter effects. Third, the resonance will only oc- 
cur for sgn(A) = -fl, whereas for sgn(A) = —1 
there will be an antiresonance with suppressed 
oscillation amplitude; c/., Eq. (1^ . This implies 
that a resonant transition occurs for neutrinos 
and Am^ > 0, or anti-neutrinos and Am^ < 0. In 
fact, one can use this effect to measure the mass 
ordering with high sensitivity, because for strong 
matter effects the rates will be strongly affected 
by the sign of Am^ . And fourth, we have learned 
above that the two-flavor probabilities should be 
CP invariant, whereas they are not in the pres- 
ence of matter, i.e., in Eq. (I?7|) . For a CP in- 
variant problem in matter, one would have to 
CP-conjugate the matter potential as well, which 
means that the Earth matter would have to be re- 
placed by antimatter (which is, of course, impos- 
sible). Therefore, matter effects violate the CP 
(and even CPT) symmetries in an extrinsic form. 
In any realistic experiment with strong matter 
effects, the CP violation has therefore to be ex- 
tracted from a convolution of the intrinsic (from 
5cp) and extrinsic (from the matter potential) 
CP violation, which implies that Eq. (jlip is not a 
very good description of CP violation in neutrino 
oscillations if matter effects are present. 

3.2. Three-flavor effects 

For the illustration of three-flavor effects, the 
most relevant oscillation channels will be the 



~^ (o"" — > t'/i) channels. In Eq. 
we have learned that Pefi — Pf_ie — for at- 
mospheric experiments to a first approximation, 
which means that deviations from zero will be 
driven by 6'i3 and the solar oscillation contribu- 
tion. In order to switch these effects on, it is 
therefore appropriate to expand these appearance 
probabilities to second order in sin 2^13 and the 
hierarchy parameter a = ATn\i/ Am\i ~ 0.03 
as |33l34l35) 



P, 



2^ sin2[(l-i)A3i] 



:^ sin^ 2^13 sin^ ^23- 

(1 - AY 

± asin20i3 sin(5cp sin2(?i2 sin2(?23 

xsin(A3r)5^^^^^*^4M 
A {l-A) 

+ a sin 2^13 cos (5cp sin 2^12 sin 26*23 

sin(iA3i) sin[(l - i)A3i] 



X cos(A 



31J 



A 



(1~A) 



a cos 023 sin 20i2 



sin^(AA3i 

i2 



(31) 



Here the sign of the second term refers to neutri- 
nos (plus) or anti-neutrinos (minus). Note that 
the sign of A, defined in Eq. (1^^ . depends on neu- 
trinos or anti-neutrinos as well. The T-inverted 
probability P^e, however, can be obtained from 
Eq. (|3ip by changing the sign of the second term 
only. 

From Eq. (PT|) . we can immediately read off 
that all of the interesting quantities ^13 , the mass 
hierarchy (by the effect in A), and (5cp (by the 
second and third terms) can be measured in prin- 
ciple if the spectral dependence of the different 
terms can be disentangled. However, because 
of the complex parameter dependence and mat- 
ter effects, continuous correlations and several 
discrete degeneracies remain in the parameter 
space even if both neutrinos and anti-neutrinos 
are used: The (5, ^13) [36], sgn(Am|j) [37], and 
(^23, 7r/2 — ^23) [3E] degeneracies, i.e., and overaU 
"eight-fold" degeneracy [S^. Using enough en- 
ergy resolution, different baselines, different os- 
cillation channels, or more statistics, the correla- 
tions and degeneracies can be resolved. One ex- 
ample is the condition sin(^A3i) = in Eq. ([5T]) . 
which makes the second to fourth terms disap- 
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Figure 4. Discovery reach for CP violation (3cr) 
for the IDS-NF neutrino factory [41] as a function 
of (true) sin^ 26*13 and (true) Scp (upper panel) 
or fraction of (true) Scp (lower panel). Normal 
mass ordering assumed. 



pear, and leads to a clean measurement of ^13 
and the mass hierarchy. This condition evaluates 
to V^Gf ne L = 2tt independent of neutrino en- 
ergy and oscillation parameters, or L ~ 7 500 km, 
the so-called magic baseline |40| . It is in general a 
good strategy to combine a shorter baseline with 
weaker matter effects in order to measure CP vio- 
lation with a longer baseline with stronger matter 
effects to measure the mass hierarchy. 

Let us take a closer look at the discovery reach 
for CP violation. A discovery of CP violation will 



be made if all CP conserving solutions Sqp = 
and TT can be excluded at a certain confidence 
level for an arbitrary (allowed) choice of the other 
oscillation parameters. In practice, one marginal- 
ized over the other parameters. From Table [Tl we 
know that the solar and atmospheric parameters 
are already very well known, whereas sin^ 2^13 < 
0.1 and dcp are unknown quantities. The per- 
formance of any future experiment will however 
depend on sin^ 2^13 and Sep within their allowed 
ranges, i.e., the values which Nature has actu- 
ally implemented. These values are often referred 
to as true values or simulated values, and cor- 
respond to data in an existing experiment. For 
the CP violation measurement, sin^ 2^13 and Sqp 
(and the mass ordering) are the critical parame- 
ters which determine the actual experiment per- 
formance. Consequently, any future experiment 
should operate in a (true) ^13 and Sep range as 
large as possible. For the quantification of the CP 
violation performance one therefore often shows 
the region in sin^ 2^13 and Sep where CP vio- 
lation will be discovered, as illustrated in Fig. |4] 
(upper panel) for a neutrino factory. Obviously, if 
sin^ 2^13 is too small, the second and third terms 
in Eq. (PT|) will be too small, and no CP vio- 
lation will be discovered. If Sep is too close to 
one of the CP conserving values, CP violation 
cannot be discovered either. A different repre- 
sentation of the CP violation discovery potential 
is shown in the lower panel of Fig. S) In this 
case, for each sin^ 26*13, the sensitive regions are 
stacked, and the fraction of Sqp for which CP 
violation will be discovered is shown. For exam- 
ple, for sin^ 26*13 = 10"^, CP violation wih be 
discovered for about 80% of all possible values 
of Sep. This representation turns out to be use- 
ful for experiment performance comparisons, see, 
e.g., Ref. [42]. Note that in Fig. E] the combi- 
nation of 4000 km and 7500 km baselines is used 
to resolve the degeneracies and to disentangle the 
intrinsic from the extrinsic (matter effect) CP vi- 
olation, which would otherwise lead to irregular- 
ities in the discovery regions. 

The mass ordering is mainly measured by the 
first term in Eq. (|3T]). which means that the sen- 
sitivity depends on sin^ 26*13. For the sake of 
simplicity, let us chose the magic baseline, where 
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only this first term survives. Then Eq. (1311) re- 
duces to the two- flavor limit in Eq. (|27l) . apart 
from the factor sin^(?23- In this case, for small 
enough 6*13, the parameter mapping in Eq. (1^51) 
is given by ^ ~ 1 — A, and the resonance con- 
dition corresponds to A 1. At the resonance, 
Pefj. oc cx L^, which compensates for the 
flux drop of the event rates. Therefore, the event 
rates stay almost constant for a wide baseline 
range; c/.. Fig. 1 in Ref. [33]. Since the effect is 
opposite for the anti-resonance (such as the other 
mass ordering), long baselines and high enough 
neutrino energies covering the matter resonance 
energy are beneficial for the discrimination of the 
mass ordering. 

3.3. Future experiment classes 

Here we consider future reactor and 
accelerator-based long baseline experiments to 
find the unknown quantities, ^13, the mass or- 
dering, and CP violation. The experiments are 
typically classified by their neutrino production 
mechanism: 

Reactor experiments with two detectors 

use neutrinos produced by beta decay in nuclear 
fission reactors, such as nuclear power plants. 
They are similar to the reactor experiments from 
the last section, measuring 9i3 in the short base- 
line limit of Eq. (j2ip . As a major improvement of 
the CHOOZ experiment, additional near detec- 
tors help to better control systematics, such as 
the normalization of the flux. Examples for these 
experiments will be Double Chooz [44 and Daya 
Bay [45]. 

Superbeams follow the technology of conven- 
tional beams producing neutrinos by (mostly) 
pion decays. Compared to the conventional 
beams, the proton beam intensity on the tar- 
get will be higher, and the detectors will be 
more massive. In addition, the off-axis tech- 
nology [46] is typically used, which means that 
the main detector is placed slightly off the main 
beam axis to reduce the beam energy and to 
over-proportionally reduce backgrounds intrinsic 
to the beam. Examples are the currently planned 
T2K ("Tokai to Kamioka") [47] and NOvA 
("NuMI Off-axis Neutrino Appearance") [35] ex- 
periments. 



Superbeam upgrades are more speculative 
ideas to push the conventional technology to its 
limits. This includes extremely high thermal 
powers in the target, and detector masses in the 
megaton class (for water Cherenkov detectors). 
One typically distinguishes two categories: nar- 
row band beams are based on the off-axis technol- 
ogy, whereas wide band beams are using a detector 
operated on the beam axis. Note that "narrow" 
and "wide" refer to the broadth of the energy 
spectrum here. There are now many ideas under 
discussion and evaluation. Examples for narrow 
band beams are upgrades for the T2K experiment 
using a megaton-size water Cherenkov detector 
(T2HK - "Tokai to Hyper-Kamiokande" ) [37]), 
or even splitting this detector mass between sites 
in Japan and Korea (T2KK - "Tokai to Kamioka 
and Korea") [49]). An example for a wide band 
beam is an on-axis beam from Fermilab (USA) 
to an Deep Underground Science and Engineer- 
ing Laboratory (DUSEL) in the Homestake mine 
(South Dakota, USA) [5D1. 

Beta beams produce neutrinos by beta de- 
cays of boosted isotopes in straight sections of a 
storage ring |51| . Compared to the superbeams, 
one has a flavor-clean electron neutrino beam 
with very predictable spectrum. However, the 
ion source has to produce enough radioactive ions 
per time frame, and a relatively large accelerator 
has to boost them to their target energies. This 
approach is currently under study from both the 
experimental and theoretical point of view, such 
as within the EURISOL ("European Isotope Sep- 
aration On-Line") design study [52]. 

Neutrino factories produce neutrinos by 
muon (anti-muon) decays in straight sections of a 
storage ring |33I53I54|55] . In this case, the spec- 
trum from the purely leptonic muon (anti-muon) 
decays is very well known, but not flavor clean: 
Since both (D^) and De (i^e) are produced si- 
multaneously in equal amounts (c/., seond line of 
Eq. ([3])), charge identification of the secondaries 
is required in the detector to distinguish the orig- 
inal flavors. In principle, the muon production 
is technologically straightforward, but the muons 
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have to be collected, coolecfl, and accelerated fast 
enough before they decay. This technology is cur- 
rently under investigation in the international de- 
sign study for the neutrino factory IDS-NF |41) . 
Very interestingly, a neutrino factory uses in part 
the same technology required by a muon collider, 
which means that it could be a first step towards 
such an experiment. The CP violation discovery 
reach of the IDS-NF neutrino factory is shown 
in Fig. m CP violation could be measured more 
than three orders of magnitude in the oscillation 
amplitude sin^ 2^13 beyond the current bound. 

3.4. Simulation of future experiments 

Typical questions regarding the optimization of 
future experiments concern the status quo at the 
time of the decision, the timescales of different 
experiments, the comparison of experiments, and 
the complementarity of the information obtained. 

As far as the optimization of individual exper- 
iment classes is concerned, one has to distinguish 
between green-field scenarios and site specific pro- 
posals. The first can, in principle, be attached to 
any major high energy laboratory (with possibly 
substantial extra effort), whereas the latter de- 
pend on a specific site, where some components 
may be already available. The study of green- 
field scenarios is especially useful to identify the 
optimal setups from the physics point of view, 
and to quantify the tradeoff for specific sites. Of 
course, the further in the future an experiment is, 
the more green-field the considered scenarios will 
be. 

As far as the quantities of interest for the opti- 
mization are concerned, there are typically: 

1. The energy of the parent particles, such as 
the 

• Ions for a beta beam (quantified by the 
boost factor 7) 

• Muons for a neutrino factory (quanti- 
fied by the muon energy E"^) 

• Pions/kaons for a conventional beam 
(typically quantified by the proton en- 
ergy Ep, where the pious and kaons 

^Here "cooling" refers to producing bunches small enough 
in phase space to be acceptable for the accelerator. 



are produced by the interactions of the 
protons with a solid target). 

2. The baseline L, i.e., the distance between 
source and detector. 

3. The integrated luminosity, which is propor- 
tional to the 

number of useful parent decays x running 
time X detector mass. 

4. Detector properties, such as efficiencies, en- 
ergy resolution, and the ability to measure 
the charge of the secondary particle. 

5. Systematical errors (and their treatment). 

6. Different parent particles used (such as dif- 
ferent isotopes for beta beams). 

7. The addition of other oscillation channels. 

8. The off-axis angle (for superbeams). 

9. Potential hybrids of different experiments. 

Whereas the baseline and parent energy can be 
easily optimized for, one can, for example, only 
optimize for systematical errors or detector prop- 
erties in this framework by identifying the quan- 
tities critical to the physics output, which is of 
interest for the experimentalists to focus their re- 
sources. 

For the simulations, often the publically avail- 
able GLoBES ("General Long Baseline Experi- 
ment Simulator") software |57I58| is used. This 
is a multi-purpose software for the simulation of 
individual long-baseline and reactor neutrino os- 
cillation experiments, as well as for the global 
analysis of multiple experiments. It includes the 
treatment of statistics, systematics, correlations, 
and degeneracies. It consists of two major compo- 
nents: Abstract Experiment Definition Language 
(AEDL) describes individual experiments using 
plain text files, and a user interface (C library) for 
the analysis, which loads one or more AEDL 
files and provides the functionality for the sta- 
tistical analysis. The separation between AEDL 
and the user interface makes GLoBES an inter- 
esting tool for both the experimentalist and the- 
orist. For example, the theorist may use pre- 
defined files for the simulation of new, potentially 
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Figure 5. CP violation discovery reach (90% CL) simulated for the T2K, NOj^A, Double Chooz, and 
Daya Bay experiments as a function of the year (left panel) and for different experiments separately (right 
panel). Inverted mass ordering assumed. Figure from Ref. [56 . 



interesting physics effects. The experimentalist, 
on the other hand, can quickly test the impact 
of modifications in the experiment definition on 
physics. Note that GLoBES is not meant to re- 
place a full Monte Carlo simulation of the experi- 
ment, but has to be understood as a tool to iden- 
tify the key parameters and critical factors for es- 
pecially future experiments. For example, the de- 
tector is simulated by an effective response func- 
tion. This response function can be used from 
Monte Carlo simulations as an input for GLoBES. 

AEDL describes an experiment, such as by 
source type and spectrum, matter density pro- 
file, cross sections, detector properties (efficien- 
cies, energy resolution, backgrounds), and sys- 
tematics. It uses three building blocks: A channel 
links a produced flavor state with a certain flux, 
via the oscillation physics, to the detection with a 
speciflc interaction type and the respective cross 
sections. It results in the event rate of this inter- 
action type. A rule combines the event rates from 
different channels, which can either be signal or 
background for that rule, with a speciflc system- 
atics; it results in a x^- An experiment contains 



one or more rules, which are combined to the to- 
tal x^. It shares certain characteristics among the 
rules, such as baseline and matter density profile, 
but not the systematical errors. For example, 
a simple neutrino factory may store , which 
leads to Ve and Df^ in the beam. A signal channel 
might be i^e — >■ '^/i, which can be combined into 
an appearance rule with the background chan- 
nel — ^ (for the charge mis-identified events) 
leading to a ■ An experiment may contain more 
such rules, such as for different appearance and 
disappearance channels and different polarities of 
the initial muons. 

One or more descriptions of experiments can 
be loaded by the C user interface. This interface 
provides the functionality to extract physical in- 
formation from the simulated event rate spectra. 
For example, it allows for sections and projec- 
tions (marginalizations) of the multi-dimensional 
fit manifold, i.e., it allows for the inclusion of cor- 
relations and degeneracies. Of course, one can 
also obtain low-level information, such as oscilla- 
tion probabilities and event rates. New features 
in GLoBES 3.0 are the fully customizable system- 
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atics interface supporting multiple sources and 
detectors (such as for reactor experiments), and 
the fully customizable external input to be added 
to the before marginalization. Heart of the 
analysis is the oscillation and rate engine, which 
includes a full three-flavor treatment, the use 
of arbitrary matter density profiles, and an ex- 
tremely high numerical efficiency with specifically 
designed numerical algorithms, such as Ref. [55] , 
The oscillation engine can be modified as well, 
which allows for the simulation of new physics 
effects. 

The results obtained with the GLoBES soft- 
ware can nowadays be found in the core of 
many long baseline experiment studies and ma- 
jor collaborative efforts, such as the international 
neutrino factory and superbeam scoping study 
(ISS) [12] and the US long baseline neutrino ex- 
periment study I50| , among many others. In ad- 
dition, application software, such as for multi- 
parameter marginalization, is available |60) . 

An example for such a simulation projecting 
the CP violation discovery reach is shown in Fig. [5] 
(for details, see Ref. [SS]). The exact sensitivity 
will depend on a number of factors, such as in- 
dividual experiment operation plans and data re- 
leases. However, such a simulation may give an 
idea of the information expected at a certain time. 
For example, one can read off from Fig. [5] that 
at the 90%CL, CP violation will be discovered 
for at most 50% of all values of Sqp until about 
2019. Given the low confidence level and the low 
parameter space coverage, a new generation of 
experiments may be necessary for a high confi- 
dence CP violation discovery. Comparing Fig. [5] 
to Fig. m (lower panel) for a neutrino factory, the 
latter could be such an instrument. 

3.5. New physics searches 

Apart from standard oscillation measurements, 
future neutrino oscillation facilities will be used 
for new physics searches. Here we show a few 
possible effects to be tested at the example of a 
neutrino factory, especially with the help of near 
detectors. 

If the new physics comes from heavy mediators 
integrated out at low energies, which applies to 
TeV- or GUT-scale physics, the effects can be de- 



scribed by a tower of effective operators using the 
SM fields |2l61l62j 

= -S^SM + -^cff^ + -^clr^ + ■ • • (32) 

with 

^,^«(x^O^ (33) 

which are invariant under the SM gauge group. 
Here A is the new physics scale. The lowest or- 
der addition to the SM is the Weinberg operator 
in Eq. ([2]), leading to Majorana neutrino masses. 
Therefore, it may be the next logical step to dis- 
cuss the implications of the higher dimensional 
operators. Assuming LHC-scale physics at A ~ 
1 TeV, one can generically estimate that the c? = 6 
operators are suppressed by (100 GeV/1 TeV)^ ~ 
0.01 and that the d = 8 operators are suppressed 
by (lOOGeV/lTeV)" ~ 0.0001 compared to the 
SM. In neutrino oscillations, typically percent 
level effects may be observable, i.e., the effects 
from d — 6 operators, whereas the effects from 
higher dimensional operators will be very hard to 
access. Therefore, we focus on d = 6 operators in 
the following which are generated at tree-level. 

The first class of effective d = 6 operators of 
interest are so-called non-standard interactions 
(NSI) 

^^^,^ = 2V2GFe;] iv^Y^LVa) (^VPl/a^) , 

(34) 

where i denote the charged leptons. Here Gf 
is the Fermi coupling constant and Pl and Pb, 
are the left- and right-handed (chiral) projection 
operators, respectively. Such operators lead to 
NSI matter effects adding to the Hamiltonian in 
Eq. ^ (for 7 = (5 = e in Eq. 

^ ee ^efi ^er \ 

ie"el.r ^Tr ■ (35) 

Note that e is the strength of the NSI effect rel- 
ative to the SM matter effect. In addition to the 
propagation in matter, the production or detec- 
tion processes can be affected by NSI. The neu- 
trino states produced in a source and observed 
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at a detector can be treated as superpositions of 
pure orthonormal flavor states j63i64i65j : 

\0 = E ^UH) ^ (36) 

a— e,/i,r 

For instance, for neutrino production by muon 
decays at a neutrino factory, one obtains e^^ for 
a = 5 = e and 7 = in Eq. ([M)) . Note that these 

and e'' are process dependent quantities. 

In writing down Eq. p4p . we do not require 
gauge invariance. If SU(2) gauge invariance is 
imposed at the effective operator level, typically 
charged Icpton flavor violating processes will be 
induced because the neutrinos come together with 
charged leptons in SU(2) doublets. If it is re- 
quired that all the charged-lepton processes van- 
ish, only the NSI operators made out of four 
lepton doublets survive which are antisymmet- 
ric in the flavor indices, i.e., a 7^ 7 and 13 ^ 5. 
Such operators can be naturally realized in the- 
ories with an SM SU(2) singlet singly charged 
scalar |66|67|68] . These models are, however, 
strongly constrained otherwise, such as by lep- 
ton universality tests [60^. Therefore, it is diffi- 
cult to find a model for large NSI from leptonic 
d = 6 operators, and for higher dimensional op- 
erators a model cannot be easily found without 
cancellations |69I70| . In summary, the model- 
independent NSI bounds are rather weak [71j 
and deserve a further test at future experiments. 
However, the prospects for c? = 6 operators (gen- 
erated at tree-level) are not very good from a 
theoretical perspective, which means that large 
NSI should come from higher dimensional or loop- 
induced operators. But these are generically ex- 
pected to be much smaller than the tree-level 
d — Q contributions, quite likely beyond the reach 
of future experiments unless the new physics scale 
is very close to the EWSB scale. 

Using neutrino factory near detectors, espe- 
cially and are interesting to be tested, 
because the neutrino factory beam does not con- 
tain tau neutrinos. The expected sensitivity for 
an OPERA-like detector is at the level of 10~^ 
to 10-"* at the 90%CL [72], maybe at the level 
where one may expect some d = 8 contributions. 



For matter NSI, the expected sensitivity for the 
NSI including the tau sector is about 10~^ at Str 
from the long baselines [T^, which is at least an 
order of magnitude beyond the current model- 
independent bounds, but maybe too large for sus- 
pecting a d — 8 contribution. Assuming d = 6 op- 
erators without charged lepton flavor violation, 
certain correlations between source and matter 
NSI are present [70], which lead to an enhanced 
sensitivity j72| . However, in this case, the sensi- 
tivity has to be compared to the model-dependent 
bounds, which it exceeds only by about a factor 
of two fT?. 

Another class of effective d — 6 operators are 
coming from integrating out heavy fermion fields, 
leading to non-unitarity (NU) of the PMNS ma- 
trix. Such fermions are often introduced in seesaw 
mechanisms at the TeV scale. In general, gauge 
invariant theories extending the SM with the tree- 
level exchange of heavy neutral fermions result in 
a dimension-six operator of the form |75I76| 

^nU' = (LcH) i^ (h^Lp) (38) 

with H = ia^H. Re-diagonalizing and re- 
normalizing the kinetic terms of the neutrinos, 
one has an effective Lagrangian 

I COS uy/ 

+ H.c. (39) 

with modified couplings to the W and Z bosons. 
Here N is an effective (non-unitary) mixing ma- 
trix which can be parameterized by 

= (!+£)[/. (40) 

Because neutrino oscillations are typically tested 
at energies below the gauge boson masses where 
the gauge bosons are effectively integrated out, 
the modified couplings effectively lead to non- 
standard four fermion interactions of the type 
Eq. ([34|) . However, the source, detector, and mat- 
ter NSI are correlated in a particular, fundamen- 
tal way (i.e., process- independent), which leads 
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Figure 6. Exclusion limit for sterile neutrinos 
measured in the two-flavor limit of Pee for sev- 
eral near detector distances d (distance to the end 
of the decays straight) and ideal, systematics-free 
near detectors at a neutrino factory (90% CL, 2 
d.o.f.; two near detectors in front of straights). 
The dashed curves illustrate the effect of includ- 
ing the averaging over the decay straight, whereas 
the solid curves are without this averaging. The 
fiducial detector masses are fixed to 200 kg. Note 
that there is no systematics included in this fig- 
ure. Taken from Ref. [78] . 



to an enhanced sensitivity. Especially the tau 
sector benefits from a near detector [77]. In fact, 
at the neutrino factory NSI from c? = 6 effective 
operators and NU lead to very similar effects for 
e^r, because the correlation between source and 
matter NSI is basically the same [74]. However, 
note that the source NSI are process-dependent, 
whereas the source NU are fundamental, which 
means that with the help of a superbeam-based 
experiment the effects could be disentangled, at 
least in principle. 

If the neutral fermion fields are light enough 
to be produced in the neutrino oscillation exper- 
iment, they enter Eq. Q directly as sterile (not 



weakly interacting) states. Although one does 
not expect a major contribution of these sterile 
states anymore, such as to describe the LSND 
anomaly [23! , small ad- mixtures of sterile neu- 
trinos are not ruled out and are a clear signa- 
ture for new physics. If the sterile neutrino mass 
splitting is significantly above Am|]^, i.e., A41 ^ 
A31 ^ A21, sterile neutrinos are best searched 
for at short baselines where A31 ^ A21 ^ 0. 
Therefore, a prominent model-independent way 
to search for sterile neutrinos is testing Eq. ([14)) 
in various oscillation channels at baselines where 
standard oscillations cannot have developed yet. 
A number of experiments have done such tests in 
the past, such as NOMAD [TH] and CHORUS [10] 
including v^- appearance. Similar tests can be per- 
formed at a neutrino factory. For example, take 
the electron neutrino disappearance probability 

Pee-l-sin2 20sin2(^^J^) . (41) 

Then one typically shows exclusion limits of the 
form in Fig. (6] For each curve, the r.h.s. is ex- 
cluded at the given confidence level. In this case, 
6* = is simulated in the data, and the Q and Am^ 
in the figure correspond to the fit values. Ob- 
viously, the optimum sensitivity (peak) depends 
on the choice of the baseline L, related to d in 
this case. Here d is the distance to the end of 
the decay straight: since the neutrino factory is 
a line source, the baseline L is not defined for 
short distances. The effect of the averaging over 
the line source is shown as dashed curves. The 
sensitivity typically breaks away for small Am^, 
because for too small arguments of the sine in 
Eq. (|4T]) the oscillation does not develop and B 
cannot be measured. Maximal sensitivity in the B 
direction is obtained at the first oscillation maxi- 
mum A = /S.m? L / [AE] ~ 7r/2, where the spectral 
effect can be easiest measured. Then higher os- 
cillation maxima are visible, until the oscillation 
averages out. For large Am^, typically the total 
event rate or systematics limits the sensitivity. 
The dependence of the sensitivity on the baseline 
is characteristic for the sterile neutrino example, 
whereas for the NSI and NU near detection the 
baseline choice only affects statistics. Note that 
such figures are not only used for new physics 
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searches, but also for the Am|]^-0i3-exclusion re- 
gion; see, e.g., Fig. 3 in Ref. [T7]. 

In summary, future facihties can be used 
to constrain a number of new physics effects, 
where we have only shown some examples here. 
Whereas the physics from higher dimensional op- 
erators may be most interesting in the context of 
LHC physics, the sterile neutrino example illus- 
trates that also the location of a near detector 
system of future experiments may be important, 
and needs to be taken into account in the opti- 
mization. 

4. Neutrino oscillations in the Sun 

Neutrino astronomy is an emerging field of neu- 
trino physics, which is so far based on the ob- 
servation of solar and supernova neutrinos. In 
2002, Raymond Davis Jr and Masatoshi Koshiba 
received the Nobel prize in physics "for pioneering 
contributions to astrophysics, in particular for the 
detection of cosmic neutrinos" . On the one hand, 
Raymond Davis Jr and his collaboration detected 
over three decades 2000 neutrinos from the Sun, 
which is an important piece of evidence for nu- 
clear fusion in the Suns interior. On the other 
hand, Masatoshi Koshiba and collaborators de- 
tected on February 23th, 1987 twelve of the 10^^ 
neutrinos which passed their detector from an ex- 
tragalactic supernova explosion. Both of these 
observations can be regarded as the foundation 
of neutrino astronomy. 

Especially neutrino oscillations in the Sun have 
also contributed to the measurement of the neu- 
trino properties. Nowadays they still provide the 
most accurate measurement of the solar mixing 
angle 9i2 . Neutrino oscillations in the Sun and in 
supernovae can, however, not be treated within 
the framework of the previous section, because 
the matter density is varying along the propaga- 
tion path and not constant. In order to illustrate 
the differences to constant matter, we focus on 
the two-flavor case in the Sun in this section, and 
show how it has lead to the measurement of 9i2 . 
The description of neutrino oscillations in super- 
novae is more complicated, because the high neu- 
trino densities lead to neutrino self interactions, 
which again lead to collective phenomena (see. 



e.g., Refs. j81|82j ). Note that supernova neutri- 
nos might also be used for the determination of 
neutrino properties, such as the observation of 
SN1987A has already lead to a bound for neu- 
trino lifetime. However, for neutrino masses and 
mixings, the predictions strongly depend on the 
parameters of the source and the availability of 
detectors. 

Neutrinos are assumed to be produced in the 
deep interior of the Sun, oscillate within the Sun 
until they reach its surface, and then propagate 
as mass eigenstates to the Earth. Therefore, 
there are no neutrino oscillations between Sun 
and Earth in our current understanding. As we 
shall see below, this can be naturally understood 
in terms of mass eigenstates emitted from the 
Sun. However, even if there was a superposi- 
tion of states emitted, coherence between Sun and 
Earth would be eventually lost, and mass eigen- 
states arrived at the Earth's surface. Therefore, 
we only deal with neutrino oscillations within the 
Sun in this section. Note, however, that the 
neutrinos start to oscillate again once they enter 
Earth matter. Therefore, if they pass substantial 
Earth matter before they are detected, i.e., they 
are detected on the "night" side of the Earth com- 
ing from below, some oscillating effect may be vis- 
ible. The difference between direct detection and 
detection after the propagation in Earth matter 
is also called day-night- effect. This effect has not 
been observed yet, which is not surprising: Ap- 
plying the solar parameters to Eq. ([30)) . we obtain 
a resonance energy of a few hundred MeV. Solar 
neutrinos only extend up to about 18 MeV, which 
is far below this resonance. Therefore, only small 
effects can be expected. In supernova neutri- 
nos, however, there is a tail of neutrinos extend- 
ing to much higher energies 0(100 MeV), where 
the Earth matter effects may be visible (see, e.g., 
Refs. [83|84j for different applications). 

Let us now first recall the differences between 
neutrino oscillations in vacuum and matter. In 
vacuum, we have - cf., Eq. (j25|) 

n - C/HdiagC/^ (42) 
1 / —Am2i COS 2^12 Am|jsin26'i2 \ 
^ Ie\ Am|i sin 2012 Ato^^ cos 26*12 /' 

where the arrow refers to the subtraction of an 
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overall phase factor. Here 
U 



cosc/12 sinyi2 
— sinf?i2 cos 6*12 

— ( ^ 

2£' V Am|i 



(43) 



Similarly - c/., Eq. (j26p ~ we have in matter 

n = UTidi^gU^ (44) 
1 / -Am2cos2^ Am2sin2^ 



with 



U = 



1 

Ye 





sin 20 


cos 6* 


sin 


— sin 


9 cos 







\ 


Atti^ 



(45) 



using the parameter mapping in Eq. In this 

case, the eigenstates of the Hamiltonian and the 
flavor eigenstates are connected similar to Eq. Q , 
with U replaced by U. Therefore, one refers to 
the eigenstates of the Hamiltonian as mass eigen- 
states in matter. In constant matter density, the 
states can be propagated using the evolution op- 
erator 



(46) 



because the Hamiltonian is not explicitely time- 
dependent. In varying matter density, however, 
the full Schrodinger equation has to be used 



dx 



(47) 



AE 



where and V'x are the amplitudes for being 
or Vx- Here Vx is a superposition of v^^ and i^r, 
because these are maximally mixed. For example, 
in the beginning, we have electron neutrinos, and 
therefore V'e = {ve\vin) = 1 and V's = {vx\vin) = 
0. As the next step, one uses a transformation to 
the mass eigenstates in matter (pi 



-ipx 



u 



(48) 



Applying this transformation to Eq. (|47l) . one 
finds 



dx 



AEi 



dx 



dx 




which is a coupled differential equation system. 
One can easily see that the differential equations 
decouple if the diagonal entries in the matrix are 
dominant, i.e.. 



IAto^I > 4:E 



de_ 

dx 



(50) 



This electron density profile-dependent condition 
(or derivations from it) is also called the adia- 
baticity condition, and neutrino oscillations in the 
Sun can be described in the adiabatic limit to a 
good approximation. Using this condition, each 
of the differential equations can be easily solved 
in order to obtain 4>i{x) = exp{i^i{x)) 4>i{0) with 
a phase factor ^i{x) depending on the matter den- 
sity profile. 

Let us now discuss the simplest case A ^ 1 
in Eq. (j29l) at the production point. That im- 
plies both high enough densities at the production 
point and large enough neutrino energies. Then 
we can read off from Eq. ((28|) that sin 20 — > or 
6 = tt/2 (which is the solution for Am2i > 0). At 
the production point a; = 0, where the neutrinos 
are produced as electron neutrinos, we therefore 
have 



V'(.T = 0) = 

coa9 \ i>i 
sinf 



1 






(51) 



At an arbitrary point x after propagation, we find 
(x) = ( \ (x) (52) 



-ipx 





pii2{x) 
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Figure 7. Electron neutrino disappearance prob- 
ability for solar neutrinos as a function of energy 
in the perfectly adiabatic limit. The thick curve 
represents the current best-fit values of the solar 
parameters, the thin curves limit the 3cr allowed 
range for sin^ 9i2 {cf., Table[l]). The dashed curve 
shows the solution for the inverted mass ordering. 



In the Sun, the electron density drops approxi- 
mately exponentially as 



.o.f.(53) 



ne{x) = ne(0) exp ( - — 



This means that for a:; — >■ oo, the electron density 
drops continuously to the vacuum density, and we 
have 



{x ^ oo) = U 







sm6'i2 
cos Oi2 



Finally, we obtain for the electron neutrino dis- 
appearance probability 



Pee — 



(1 o)( 



(t oo) 



(55) 

i.e., the phase factor, which we have not com- 
puted explicitely, drops out. 

The factor sin^ 9i2 describes the disappearance 
of electron neutrinos from the Sun through neu- 
trino oscillations in the Sun's interior at high 
enough energies, such as measured by the SNO 
experiment. In practice, the condition A ^ 
1 used for this derivation only holds for E > 
10 MeV, whereas for very small energies E < 
1 McV one basically obtains the vacuum result 
{A ~ 0) with averaged oscillations {cf., Eq. (fT2|) 
applied to Eq. ([T4|) ) 



1 



1 



■ sin-^ 26*12 . 



(56) 



For the intermediate case, one finds a continu- 
ous transition depending on the size of A. For 
Atoji > 0, the effective mixing angle 9 in Eq. (|28p 
starts at the vacuum value 9 at low energies, in- 
creases to 7r/4 at the resonance energy (where 
sin 20 — 1), fiips the octant and increases fur- 
ther to tt/2. For Ato2i < 0, there is no reso- 
nance, and the angle decreases continuously from 
the vacuum angle 9 to 0. In this case, sin 29 
at high energies evaluates to 9 ^ 0. One ob- 
tains (j){x = 0) = (1,0)"^ in Eq. leading to 
Pee = cos^ 9i2. We illustrate the transition prob- 
ability as a function of energy in Fig. [7] for the 
perfectly adiabatic case and the parameters from 
Table [1] (thick curve). In this figure, the lower en- 
ergy limit corresponds to Eq. (1551) . the upper en- 
ergy limit end to Eq. (j55p . The thin curves limit 
the 3(7 allowed range for sin^ 6*12 {cf., Tabled]). 
Note that there is little dependence on Am|i in 
its 3(7 allowed range, because Am^i is very well 
measured by the KamLAND experiment. The 
dashed curve shows the Am|i < case. The en- 
ergy dependence in Fig. [7] has been measured by 
early Gallium and the Homestake experiments at 
the very left end, and by SNO at the very right 
end. Therefore, the Arn^i < case has been ex- 
cluded from the solar neutrino observations H In 



= sm U12 



"Strictly speaking, this discussion can only be done to- 
gether with the choice of the octant of 612 , since there is an 
ambiguity sin^ 612 — > cos^0i2, Am|-^ — > — Am-jj in these 
derivations, i.e., instead of changing the sign of Am|]^, one 
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the future, the BOREXINO experiment (5S| has 
the potential to improve the information in the 
intermediate energy range, and verify the energy 
dependence which is so characteristic of the so- 
lar flavor transitions. Note that the flavor transi- 
tions described in this section are also often called 
the MSW effect, named after Mikheyev, Smirnov, 
and Wolfenstein [29 30 31) . The case of constant 
matter density in the previous section is a special 
case of the general MSW effect. 

5. Neutrinos from cosmic accelerators 

High energetic neutrinos are especially pro- 
duced in the Earth's atmosphere or in man-made 
terrestrial experiments. However, above the TeV- 
boundary, also cosmic accelerators may produce 
neutrinos, see Refs. |86l87l88j for review arti- 
cles. In particular, the observation of high en- 
ergetic cosmic rays, which are believed to come 
from extragalactic sources at very high ener- 
gies, motivates the existence of such accelerators. 
If, however, these hadrons interact with other 
hadrons or photons in the source, a significant 
neutrino flux will be produced as well. Exper- 
iments such as ANTARES [8^ in the Mediter- 
ranean or IceCube [20] at the South Pole are built 
for the observation of such fluxes, they are often 
called neutrino telescopes. Known candidates for 
extragalactic accelerators as potential neutrino 
sources include active galactic nuclei (AGNs) [91] 
192193) and gamma ray bursts (GRBs) [M], see 
also Ref. [95] for theoretical considerations. If 
these sources accelerate enough hadrons, neutrino 
fluxes are to be expected, which should be observ- 
able in the neutrino telescopes. However, such 
cosmic neutrinos have not been observed yet, a 
fact which may be not so surprising from the cur- 
rent point of view. If one relates the possible 
neutrino flux from such cosmic accelerators to the 
cosmic ray flux, one obtains an upper bound on 
the neutrino flux, the so-called Waxman-Bahcall 
bound [SS]- This bound assumes that the neu- 
trons, produced in photohadronic processes, es- 
cape from the source before they decay. A differ- 

can also change the octant of 812- In cither case, there are 
two qualitatively different cases (resonance/no resonance), 
which can be distinguished. 



ent version is the Mannheini-Protheroe-Rachen 
bound [97| , which includes sources optically thick 
to neutrons (the neutrons interact before they can 
escape), and uses gamma rays as an additional 
information source. The IceCube experiment will 
exceed these bounds in the coming years, which 
means that the detection of cosmic neutrinos in 
the near future might be quite plausible. 

The observation of extragalactic neutrino 
fluxes may be interesting for different reasons. 
First of all, the fluxes are evidence for the hadron 
content in the source. Second, neutrinos di- 
rectly point to the source, unlike the cosmic rays, 
which are affected by magnetic fields, and pho- 
tons, which are easily absorbed or scattered. And 
third, such neutrino fluxes may be used for the de- 
termination of new physics properties. For exam- 
ple, neutrino properties such as neutrino lifetime 
will be tested. In this section, we will argue from 
the source via propagation to detection, focusing 
on the determination of neutrino properties via 
propagation effects. 

5.1. Neutrino production in the source 

Neutrinos are typically assumed to be pro- 
duced by pp interactions, or p7 photohadronic 
processes. These interactions lead to charged pi- 
ous, among other particles, which decay into neu- 
trinos through Eq. ([3|). The interacting protons 
are assumed to be accelerated in relativistic jets 
by Fermi shock acceleration, which leads to a 
power law proton injection spectrum. The target 
protons in pp interactions are typically introduced 
using external material, such as dust, hit by the 
relativistic outflow, which needs to be described 
by additional parameters. In self-consistent mod- 
els, photohadronic interactions with target pho- 
tons are used, which originate from synchrotron 
radiation of electrons or positrons co-accelerated 
with the protons. The synchrotron photon field 
can, to a first approximation, be computed from 
the model parameters, such as the magnetic 
field and the spectral index of the electrons 
(positrons), see, e.g., Ref. [101] for such an ap- 
proach. Alternatively, a thermal target photon 
field, such as from an accretion disk, may be used 
as a target, which again introduces new parame- 
ters. 
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Figure 8. The total P7 photo-meson cross section 
as a function of the photon energy in the pro- 
ton rest frame analog to [55] (l^barn = 10 
cm^; data, shown as dots, from Ref. [99])). The 
contributions of baryon resonances (red, dashed), 
the direct channel (green, dotted) and multi-pion 
production (brown) are shown separately. Figure 
taken from Ref. [TOD] . 
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The photohadronic neutrino production is of- 
ten described by the A-resonance approximation 
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P- 



1/3 of all cases 
2/3 of all cases 



(57) 



This description, is however, not sufhcient for our 
purposes. First of all, other processes contribute, 
as shown in Fig. [8l such as direct (t-channel) pro- 
duction or higher resonances. These other pro- 
cesses have different characteristics, such as dif- 
ferent energies where the pious are found as a 
function of the initial proton energy (different in- 
elasticities) , and different multiplicities of the pi- 
ous. For example, 7r~ are produced in higher res- 
onances in addition to the 7r+ and tt'^ in Eq. ([57]) . 
These tt" affect the neutrino-antineutrino ratio, 
which may be used to test the difference between 
pp and pj neutrino production (in the pp produc- 
tion TT"*" and 7r~ are produced in equal amounts) . 
The Glashow resonance De + e~ — > ... at 
about 6.3 PeV [10211031104] is a detection pro- 



cess suitable for neutrino-antineutrino discrim- 
ination. In addition, the high energy tail of 
the multi-pion production in Fig. [8] leads to a 
change of the spectral shape, as we will see below. 
More refined descriptions of the neutrino fluxes 
from photohadronic interactions therefore use the 
Monte-Carlo simulation of the SOPHIA soft- 
ware ^98|il05j . For time-consuming source simula- 
tions, an effective description such as Ref. |106| is 
an efficient alternative, or a simplified interaction 
model, such as Ref. |100j if the intermediate par- 
ticles, such as muons, are to be treated explicitely 
to include cooling effects. 

The individual contributions of the different 
processes to the 7r+ and 7r° spectra in a GRB 
are shown in Fig. [9] in the left and right pan- 
els, respectively. In these spectra, the cutoff at 
high energies is introduced by hand, since cool- 
ing processes are not included. The spectrum 
from resonant production, shown as thin solid 
curves, resembles the typical GRB spectra often 
shown in the literature. However, the slope of 
the plateau is increased by the fact that the cross 
section in Fig. [8] is non-vanishing for large en- 
ergies, as it can be seen by the dashed curves 
representing the multi-pion contribution, and in 
fact also depends on the extrapolation to higher 
energies. For the 7r+ spectrum, the resonant pro- 
duction only dominates in a very narrow region 
around the first spectral break, whereas the 7r° 
spectrum is dominated by the resonant produc- 
tion, mostly Eq. (|57p . The reason is that there are 
more charged pions produced in the direct pro- 
duction, which compensates for the lower cross 
sections. Therefore, if Eq. (|57|) is used to estimate 
the neutrino flux from the photon flux coming 
from the tt*^ decays, the neutrino fluxes are typ- 
ically underestimated. From the center of mass 
energy dependence of the interactions, it can be 
shown that, independent of the input spectra, the 
ratio between charged pions and neutral pions is 
not 1:2, as in Eq. dSTj), but at least 1.2:1 [TUD] . 

The neutrino production finally is given by the 
weak decays of pions and muons, as described in 
Ref. [107j . Very importantly, the muon decays 
are helicity dependent, which means that the spin 
state of the muon has to be taken into account, 
and the four muon species /zj, ^Jj, and /i^ 
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Figure 9. Contributions of resonant (thin solid), direct (dotted) and multi-pion (dashed) production for 
tt"*" (left) and 7r° (right) spectra for proton-photon interactions using the GRB benchmark from Ref. [107) . 
Figure taken from Ref. |100) . using the description of processes in there. 



have to be treated separately. 

5.2. Neutrino fluxes 

For the observed neutrino fluxes, one typically 
distinguished three different conceptual cases. A 
single source or point source flux is related to a 
particular source, such as a GRB or AGN. Un- 
fortunately, the statistics to be expected from 
a single (extragalactic) event is rather moder- 
ate, a few detected neutrinos at most; see, e.g., 
Ref. [108) . Therefore, different techniques have 
been proposed in the literature to increase the 
statistics. 

The first approach is a stacking analysis; see, 
e.g., Ref. |109| for an example. In this case, the 
expected event rates from different, but similar 
sources are added. For example, the information 
from other messengers, such as gamma rays ob- 
served by BATSE or Fermi-LAT, can be used to 
compute the predicted neutrino flux under cer- 
tain assumptions. By stacking the very few ex- 
pected neutrino events from each source, reason- 
able statistics might be obtained with a good sig- 
nal to background ratio, because time and di- 
rectional information can be used to reduce the 
atmospheric neutrino background. One problem 



of the approach in Ref. |109| is that the redshift 
of each source is needed to reconstruct the neu- 
trino spectrum from an observed photon spec- 
trum, which is only measured in a few cases. In 
addition, a few transient events seem to dominate 
the obtained neutrino prediction, which means 
that there may be selection effects. 

The second possibility to increase statistics is 
measuring the diffuse flux from all sources in the 
sky. A generic formula for the observed flux is 
given by [57] 

-{£,z) 



z C 



1 



'^ATTdLiz^- ^^^^ 

Here n is the source distribution function as a 
function of luminosity C and redshift z, d'^u/dEy 
is the single source flux, and d^ is the lumi- 
nosity distance. The energy at the source 
is redshifted by the cosmic expansion to E'^ = 
E^/{1 + z). Note that the source flux is already 
given in the observer's frame in this case. Dif- 
fuse fluxes provide the highest statistics, since all 
possible sources contribute. However, due to the 
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lack of directional and timing information, the 
signal to background ratio may be poor. Espe- 
cially at lower energies, the atmospheric neutrino 
background prohibits an extraction of the diffuse 
flux. 

In summary, a point source flux provides the 
cleanest information, because there are no av- 
eraging effects involved. However, the statistics 
is poor. Stacked and diffuse fluxes increase the 
statistics significantly. It is however unclear what 
the impact of averaging is on the discussed mea- 
surements of neutrino properties. 

5.3. Flavor composition and propagation 

The discussion of particle physics properties of 
neutrinos often is based on the flavor composition 
at the source, which may be changed by propa- 
gation effects. Generically, three different source 
classes are distinguished, where neutrinos and an- 
tineutrinos are not discriminated: 

Pion beam sources produce neutrinos of the 
flavors Ue'-v^WT in the flavor ratio 1:2:0, as 
it is expected from Eq. ([3]). 

Muon damped sources assume that the 
muons loose energy before they decay, 
which means that at high energies Vf,:v^:vT 
come in the flavor ratio 0:1:0 according to 
Eq. ®. 

Neutron beam sources assume neutrino pro- 
duction by neutron decays, which may come 
from the photo-dissociation of heavy nuclei. 
Therefore, one has Vf,:v^:Vr in the flavor ra- 
tio 1:0:0. 

Typically these flavor ratios are discussed in an 
energy-independent way. However, in practice, 
the flavor ratios change as a function of energy. 
For example, assume that neutrinos are produced 
by pion decays. In this case, the muons loose en- 
ergy by synchrotron radiation already at lower 
energies than the pions because of the smaller 
mass. This means that the pion beam source 
changes into a muon damped source at the energy 
where the muon decay and cooling timescales are 
equal [1101, which is also sometimes called muon 
damping. The actual energy dependence of the 
flavor ratios is non-trivial and depends not only 



on cooling and decay timescales, but also on the 
participating particle species and their interac- 
tions, see, e.g., Refs. |111|107I112] . On the other 
hand, the energy dependence contains non-trivial 
information on the astrophysical properties. 

The propagation of the neutrinos between 
source and detector consists of two parts. First 
of all, flavor mixing 

3 

Pa/3-^|t/„.nf/ftP, (59) 
1=1 

corresponding to neutrino oscillations with the os- 
cillating part averaged out, changes the flavor ra- 
tios. The averaging of the neutrino oscillations 
can, for example, be justified by the size of the 
production region or decoherence on the way to 
the Earth. For example, if the original neutrinos 
are produced in the flavor ratio Vf.:v^:Vr of 1:2:0, 
one will, depending on the mixing parameters, 
roughly have 1:1:1 at the Earth; see Ref. [113] 
and references therein. Note that Eq. (ICTl) is in- 
sensitive to CP violation, but it depends on the 
CP conserving part cos^cp, which may be used 
to extract some information on the phase. The 
second part of propagation effects concerns with 
neutrinos passing the Earth's interior before be- 
ing detected. At high energies, the Earth becomes 
opaque to and because of the interaction 
cross section increasing with energy. On the other 
hand, the Vt are partly regenerated in the prop- 
agation. 

5.4. Neutrino detection and flavor ratios 

Neutrinos are detected with neutrino tele- 
scopes, which are primarily designed for the 
tracking of muons in water or ice coming from 
muon neutrino interactions. The muons and sec- 
ondaries can be seen through Cherenkov radia- 
tion in photodetectors. The threshold for this 
process depends on the spacing of the photode- 
tectors. As a peculiarity, the muons do not have 
to be produced in the actual detector volume or 
even pass it to be detected, only their light has 
to be received. The detection of muons is there- 
fore described by an effective area depending to 
the actual detector geometry and the muon range, 
which is a function of energy. 
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Measuring flavor ratios of astropliysical neutri- 
nos, on the otfier liand, is much more sophis- 
ticated (see Ref. [114] for an overview). Elec- 
tron and tau neutrinos produce (apart from neu- 
tral current events) electrons and tauons, respec- 
tively, which result in showers (cascades) of par- 
ticles. These can only be identified if the in- 
teraction vertex is within the detector volume. 
Therefore, these events arc fiducial volume dom- 
inated, and the expected event rates are much 
lower. The discrimination between electromag- 
netic showers (from electrons) and hadronic show- 
ers (from tauons) is very difficult. At higher en- 
ergies, however, the tauon may live long enough 
such that its track is separable, which means 
that the tau event may be identified as by so- 
called double bang or lollipop events |102| . de- 
pending on the visibility of the first interac- 
tion vertex. In addition, the above mentioned 
Glashow resonance De + e~ W~ —j' ... at about 
6.3 PeV I102I103I104] may be used for (electron) 
antineutrino identification, whereas all other pro- 
cesses cannot distinguish between neutrinos and 
antineutrinos. 

For a simplified (energy-independent) discus- 
sion, it turns out to be useful to define observables 
which fulfill the following two requirements: 

1. They take into account the unknown fiux 
normalization. 

2. They take into account the detector prop- 
erties. 

For requirement 1), flavor ratios, i.e., the ratios 
of fluxes at the detector between different flavors, 
are the most popular choice, since in these ratios 
the flux normalization cancels. For requirement 
2), it is convenient to define with increasing level 
of difficulty for the detection the following sim- 
plified ratios: 



R 



ADet 



ADet 



(60) 



where c^J^*^* is the neutrino flux (accumulated over 
a certain energy range) of fiavor at the detec- 
tor, neutrinos and antineutrinos summed over. 
The ratio R corresponds to the ratio of muon 



tracks to showers, where electromagnetic and 
hadronic showers do not need to be discriminated 
(here the background from neutral currents is ig- 
nored). Naturally, it is fiducial volume limited, 
i.e., the number of observed showers limits the 
statistics. As a second observable, we define 



iDct 



ADct 



(61) 



which basically corresponds to the number of 
electromagnetic (from i>e) to hadronic (from Vr) 
showers, or electromagnetic showers (from Ve) to 
double bang or lollipop events (from 1^^-) if the 
T track can be resolved. Obviously, it is more 
difficult to measure S than R, such as because 
of the limited statistics of the events which can 
be uniquely identified (for example double bang 
events). Furthermore, 



T = 



iDet 



/jDet 



(62) 



represents the ratio of Glashow resonance elec- 
tron antineutrino events (described by 



j,Det\ 



to 



the muon tracks. Note that these ratios do not 
take into account neutral current processes, dif- 
ferent statistics between numerator and denom- 
inator, and backgrounds (such as from different 
r decay modes), which means that results from 
these should be interpreted with care. 

These flavor ratios have been extensively used 
in the literature to identify a number of possible 
standard and new physics effects, see Ref. |113] 
for a review. Here we focus on two applica- 
tions: Flavor mixing, which is a standard effect 
described by Eq. (|59p . and neutrino decay, which 
is too slow in the SM and therefore requires some 
new physics contribution. 

5.5. Constraints on neutrino flavor mixing 

The measurement of flavor mixings in Eq. ()59p 
has been extensively studied in the literature, see, 
e.g., Refs. pi)4'115'116'117"118"119'120'121'122J 
I123ll24ll25n"2 6 127 128 129 1 30 131 132 133.134] 
in this context. Although the uncertainties in the 
parameters are probably too large for realistic 
applications, at least for pion beam sources |121) . 
there could be interesting information in special 
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cases, such as using different sources or in com- 
bination with new physics scenarios, see, e.g., 

Ref. [ns]. 

Especially interesting may be the dependence 
of i?, which is the flavor ratio possibly easiest 
to observe, on 6*13 and (Sep |117| . For example, 
one may expand R for the different astrophysical 
sources to first order in 6'i3 |119j 

^Neutron beam _ q 36 + 0.30 ^13 COS <5cP , 
^Muon damped _ q gg „ q 53 6I13 COS (5cP , 
^Pionbeam _ q 5Q _ ^ 

(63) 

Higher order terms in are relatively small (but 
not negligible). Now compare this to superbeams, 
which are typically optimized for the the first os- 
cillation maximum. At the first oscillation maxi- 
mum and in vacuum, one finds from Eq. pt |) 

P^e - 2 ± 0.09 013 sin ,5cp , (64) 

where the plus is for antineutrinos and the mi- 
nus for neutrinos. Obviously, these experiments 
are optimized for the observation of CP viola- 
tion, which is proportional to sin 5cp (the imag- 
inary part of exp(i(5cp))i whereas the astrophys- 
ical sources can only capture the CP conserving 
part because of the averaging of the oscillations. 
Therefore, there is obviously complementary in- 
formation. 

From Eq. (|63p . we can read off that the mod- 
ulation with (5cp is smallest for the pion beam 
source, which may be the most common case. 
Therefore, one can probably not expect any in- 
formation from such a source. In the following, 
we will therefore focus on the other two cases. 
For example, consider Fig. [TUl where the com- 
plementarity among superbeams, reactor exper- 
iments, and astrophysical sources is illustrated 
at the example of a muon damped source. In 
this case, the simulated rates and simulated R 
are computed for sin^ 2^13 = 0.1 and (5cp = tt. 
The curves show the points in the sin^ 20i3-(5cp- 
plane which produce exactly the same observ- 
ables, i.e., are degenerate. Obviously, there is 
no information on <5cp from the reactor experi- 
ment alone, because the probabilities do not de- 
pend on (5cp; c/., Eq. (|2ip . Even if neutrinos and 



2n 







Superbeam, v^^^ 


X ^Superbeam, y 

N 

\ 

1 
/ 



Muon damped flux 


✓ 

\ t 





0.06 0.08 0.1 0.12 0.14 



Figure 10. Synergy among superbeams, reactor 
experiments, and astrophysical fluxes (at the ex- 
ample of a muon damped source) in the sin^ 2^13- 
^cp-plane. Shown are the curves for constant 
total rates (superbeam, reactor experiment) and 
constant R (astrophysical flux) going through the 
best-fit point sin^ 2^13 =0.1 and (5cp = tt. Figure 
taken from Ref. |119j . 



antineutrinos are used, there a degenerate solu- 
tion remains for the superbeams. The astrophys- 
ical source alone, on the other hand, cannot limit 
sin^ 26*13, because there is a strong degeneracy 
with cos(5cp; c/., Eq. (|63l) . The combination of 
superbeam or reactor experiment and astrophys- 
ical source, however, can in this case uniquely de- 
termine the point in the sin^ 20i3-(5cp-plane. This 
synergy cannot only be used for the measurement 
of (5cp, but also for the mass hierarchy, as illus- 
trated in Ref. [119;, because the degenerate solu- 
tion in sgn(Am|j^) moves in the Jqp direction as 
a function of sin^ 26*13 | 135| . Note, however, that 
this complementarity is only relevant for large 
values of sin^ 2^13 and as long as there is no infor- 
mation on cos (Sep from terrestrial sources, which 
can be obtained by wide band beams (beams with 
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Figure 11. Fit region as function of sin 29i3 
and ^cp for R, measured from a neutrino beam 
source at the level of 10%, combined with Dou- 
ble Chooz. The simulated values are chosen as 
marked by the dot. The contours are shown 
for the 1(7 (black curves, dark regions) and 90% 
(gray curves, light regions) confidence level (1 
d.o.f.). Dashed curves represent the results when 
the other (not shown) oscillation parameters are 
fixed, i.e., not marginalized over. The arrows 
mark the ranges in Sqp which can be excluded 
at the 90% confidence level. Figure taken from 
Ref. [112]. 



a broad energy spectrum), such as on-axis super- 
beams or neutrino factories. 

The complementarity is shown at the exam- 
ple of a reactor experiment, sensitive to sin^ 20i3, 
and a neutron beam source, with R measured 
at the level of 10%, in Fig. [TTJ This example 
is, of course, somewhat optimistic, since a 10% 
measurement of R requires 0(100) events already 
from the statistics point of view. However, it can 
be seen from this figure, that some information 
on Sep is obtained. Although there is no CP vio- 
lation measurement in this case because the simu- 



lated value Scp = TT conserves the CP symmetry, 
a small range of i5cp around can be excluded at 
the 90% confidence level. 

5.6. Constraints on neutrino decay 

Apart from "conventional" parameters, such as 
neutrino masses and mixing angles, a very promi- 
nent example new physics properties which can 
be tested by fiavor ratios is the neutrino lifetime; 
see, e.g., Refs. |107ll37ll36Tra5| . Phenomenolog- 
ically, from the observation of neutrinos from su- 
pernova 1987A, we know that at least one neu- 
trino mass eigenstate must be stable over inter- 
galactic distances. More stringent and explicit 
bounds can be derived from different observa- 
tions when specific decay models are assumed 
(see, e.g., Refs. |113ll39ll40j for an overview). 
For example, solar neutrinos strongly limit the 
possibility of radiative decays |141| , while for Ma- 
joron decays [1421143] explicit bounds can be ob- 
tained from neutrinoless double-beta decay and 
supernovae |144| . Purely phenomenological (i.e., 
mo del- independent) bounds are, however, much 
weaker, leaving enough parameter space for the 
decay of any mass eigenstate over extragalac- 
tic distances [145 146 147 1^. Note that these 
bounds strongly depend on the mass eigenstate 
considered. In view of the rather weak direct 
neutrino lifetime limits and the recently proposed 
unparticle models, which may lead to new mecha- 
nisms of neutrino decay [149I150I151I152) , it may 
be well motivated to study the most general phe- 
nomenological case. 

Conceptually, one may distinguish two types of 
neutrino decay |153j : Invisible decays often refers 
to the decay into invisible states, such as ster- 
ile neutrinos hardly mixing with the active ones, 
or unparticles. Visible decays refers to decays in 
observable states, such as other active neutrinos. 
For example, in a degenerate mass scheme, invis- 
ible decays may be the most plausible assump- 
tion, because the decay kinematics depends on 
the mass differences. In a normal hierarchical 
spectrum, the heavier mass eigenstate TO3 may 
also decay into one of the two lighter ones. An- 
other conceptual classification refers to the life- 
time of the neutrinos. Depending on the lifetime 
and energy (see below), the neutrinos may have 
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Figure 12. Allowed regions at 99% CL in the {R, S) plane corresponding to different complete decay 
scenarios, for a muon damped source (left panels) and a pion beam source (right panels). We assume a 
normal hierarchy. The upper panels correspond to the analysis of present data reported in Ref. [lOj . The 
other panels show the impact of 3 years of Double Chooz data taking (1.5 with near detector), assuming 
no signal (sin^ 26i3 — 0, middle panels) or a large signal (sin^ 26*13 = 0.1, lower panels). The branching 
ratio parameters have been varied as well, where applicable. Figure taken from Ref. [136) . 
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completely decayed over extragalactic distances, 
which is often called complete decays. On the 
other hand, incomplete decays show the charac- 
teristic energy dependence of decays. 

Including incomplete invisible decays, the 
propagation effects in Eq. (|59l) change into 

3 

Pa/j-5]|C/a.|'|t^ftl'e-"'*, (65) 

i=l 

where = rai/r^ is the decay parameter as a 
function of the rest frame lifetime . The char- 
acteristic energy dependence in the last term of 
Eq. (p5)) comes from the Lorentz boost of the rest 
frame lifetime — — T^E/irii = E/ai. 

For a recent discussion of incomplete decays as 
a function of energy in astrophysical fluxes, see 
Ref. |138j . If the decay of mass eigenstate j is 
complete, the exponential in the sum of Eq. (j65|) 
suppresses the jth term. For a more refined dis- 
cussion of possible combinations of neutrino os- 
cillations and decay, including visible incomplete 
decays, see Refs. [153|154j . 

For complete decays, no matter if visible or in- 
visible, it can be shown that only 2'^ = 8 pos- 
sibilities exist, because each of the three active 
mass eigenstates can be only stable or unstable. 
If all mass eigenstates are stable, we have the 
standard flavor mixing scenario in Eq. (j59p . if 
all mass eigenstates are unstable, no signal will 
be observed. Note that decay chains, such as 
3 — 2 — 1, can be integrated out in this ap- 
proach, since eventually all and 1^2 will end 
up as vi. The eight possibilities discussed here 
are shown in the plot legend of Fig. [TH on the 
right hand side. In this figure from Ref. [136) , 
a filled disk represents a stable mass eigenstate, 
and an unfilled disk an unstable mass eigenstate. 
The different regions in the panels show the al- 
lowed regions at 99% CL in the (i?, S) plane cor- 
responding to the different complete decay sce- 
narios in the legend, for a muon damped source 
(left panels) and a pion beam source (right pan- 
els). A normal hierarchy is assumed herel^l The 
upper panels correspond to the analysis of present 

^In fact, there is some dependence on the mass hierarchy 
in such figures. Although the flavor mixing is insensitive 
to the mass hierarchy, the possible decay channels depend 



data reported in Ref. [TU] . The other panels show 
the impact of 3 years of Double Chooz data tak- 
ing (1.5 with near detector), assuming no signal 
(sin'^ 26*13 = 0, middle panels) or a large signal 
(sin^ 26'i3 — 0.1, lower panels). The branching 
ratio parameters have been varied as well, where 
applicable. 

A particular measurement of R and S will lead 
to one specific point in the planes of Fig. [T^l It is 
clear that some scenarios could not be disentan- 
gled easily, even if both R and S were measured to 
high precisions. However, if, for instance, a large 
value of R is measured, then scenario #5 can be 
easily identified even without measurement of S. 
In fact, it turns out that this observation is even 
robust for arbitrary sources without Vr produced 
at the source. Therefore, the result at the end 
depends on the scenario actually implemented by 
nature. Furthermore, a future measurement of 
013 will somewhat reduce the regions, especially if 
^13 is found to be small, as it can be read off from 
the middle row of the figure. Finally, it turns out 
that for this particular application the pion beam 
source (right column) is better suited, because it 
leads to smaller regions with respect to the os- 
cillation parameter ranges. The reason is exactly 
the same as in the flavor mixing case: Accord- 
ing to Eq. (|63p . the pion beam is least affected 
by the unknown oscillation parameters (which is 
good for this application, but not so good for the 
flavor mixing measurement). 

6. Determination of the absolute neutrino 
mass scale and the nature of neutrinos 

From Fig. [21 we can read off that the spectrum 
of the neutrinos is fixed for known mass squared 
differences Am^i and Am|i, apart from an over- 
all scale, the absolute neutrino mass scale. This 
scale is determined by any of the three neutrino 
masses, such as the lightest neutrino mass mi 
(normal ordering - NO) or TO3 (inverted order- 
ing - 10). In practice, it turns out that specific 
processes or classes of experiments are sensitive 
to particular combinations of the neutrino masses 
including the mixing angles and, for the case of 

on the mass hierarchy. Therefore, the inverted hierarchy 
scenarios look somewhat different in general. 
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decay, also the Majorana phases. We em- 
phasize the relationship between the observables 
and the absolute neutrino mass scale in this sec- 
tion. For more details, see, e.g., Ref. pTj. 

The most direct access to the absolute neutrino 
mass is the test of the kinematical effect in nuclear 
beta decay, typically of tritium: 



He 



(66) 



In this case, the endpoint of the electron spectrum 
depends on the mass of the neutrinos, no matter 
if Dirac or Majorana particles, but the effect is 
tiny. The combination of neutrino masses, this 
experiment class is sensitive to, is given by 
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(67) 



which comes from the incoherent sum over the 
contributions from different mass eigenstates in 
Eq. ([TO)) . In terms of the neutrino oscillation ob- 
servables and the lightest neutrino mass mi (nor- 
mal ordering) or (inverted ordering), Eq. (1671) 
can be expressed as 
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|Am2i|c?3 + Am2,s2^c?3 (10) 

(68) 



sinOij and c,y 



cos 9ij. In degen- 
erate case, the mass squared splittings are neg- 
ligible, and the mixing effects can be neglected. 
In the normal hierarchical case (toi ~ 0), the 
relative contribution of TO3 depends on ^13. For 
sin^ 26*13 ^ 0.04, the TO3 contribution dominates, 
for sin^ 2^13 < 0.04, the TO2 contribution. For the 
inverted hierarchical case (m3 ~ 0), both mi and 
TO2 contribute almost equally, with the ratio given 
by tan^6'i2. Note, however, that the hierarchical 
cases will not be accessible by the next generation 
of experiments. The most stringent bound on mp 
was obtained by the Mainz [155] and Troitsk |156) 
experiments, about 2eV (95% CL). This limit is 
expected to be exceeded by one order of magni- 
tude by the KATRIN experiment [TFT] . 

Another very important test of the absolute 
neutrino mass scale is the Oi//3/3 decay. In this 
case, a nucleus, for which the single beta decay 



is disfavored because of the pairing interaction in 
the nucleus, decays (for example) via 



N{A, Z) N{A, Z + 2) + 2 



(69) 



without emitting neutrinos. This processes can 
be associated with the Majorana neutrino mass. 
Simply speaking, there it corresponds to two si- 
multaneous beta decay processes for which the 
outgoing neutrino of the first process becomes the 
incoming neutrino of the second processes^ This 
connection requires that the neutrino be its own 
antiparticle, i.e., has Majorana character. Be- 
cause 0iyf3/3 decay does not work in the SM, in 
which Eq. ([2]) is not present, an observation would 
probably one of the most exciting discoveries in 
neutrino physics. Note that the signature oi0v/3(^ 
decay is very different from 2i>f3/3 (two simultane- 
ous beta decays with neutrino emission). There 
are only two particles emitted in 0i//3/3, namely 
the electrons, leading to fixed sum of the elec- 
tron energies. For 2z^/3/3, one expects a continu- 
ous distribution of the sum of the electron ener- 
gies, because the neutrinos can carry away energy. 
The combination of neutrino masses measured in 
Oi^/3/3 decay is |mee|^. With the definition of the 
Majorana phases in Eq. ([TBI) it is given by (see, 
e.g., Ref. 162 and references therein) 
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TOi \Uel\ = TO1C12C13 , 

= TO2 |J7e2p = TO2 §12 Ci3 , (72) 
2 

13 ■ 



= TO3 |C/e3p = 7713 



In terms of the neutrino oscillation observables 
and the lightest neutrino mass 7771 (normal order- 
ing) or 7713 (inverted ordering), Eq. (j67l) can be 



I'^Strictly speaking, Eq. l(2j leads to Ou/Sfi decay of this 
kind. However, other effective operators may induce Ou^p 
decay, ft can be shown that whatever operator causes 
Ouf5/3 decay, also leads to Eq. ((2]l, and therefore Majorana 
neutrino masses, at the loop level; see generic argument 
in Ref. [T58l. 
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General theory prediction 



Disfavored by Ov/S/S decay ^^^M / 


Am3i<0(IO) 


>. 

D) 
ri 






mol( 




(0 




o 
o 




>. 




XI 




T3 







Arrigi >0 (NO) / 


avo 




CO 




Q 



10"" 10"^ 10"^ 10"^ 10° 



m [eV] 

Figure 13. General theory allowed region for 
|mee| in 01^/3/3 decay as a function of the lightest 
neutrino mass m.We show the currently allowed 
regions for |mee | , where the mixing angles are var- 
ied in the current 3(7 ranges from Table [T] (com- 
puted with the formulas given in Ref . [159] ) . Note 
that fixing 6*13 would result in the appearance of 
the "chimney". The limit from cosmology is as- 
sumed do be 1 eV for the sum of the three neu- 
trino masses, and the the limit from Oi^/3/3-decay 
is obtained by the Heidelberg-Moscow collabora- 
tion [160) . We fix the mass squared differences to 
their best-fit values in Table [TJ Figure updated 
from Ref. [161] . 



expressed as 

|TOee| = '7iiC22Ci3 \J ^1 + ATO21 ^12^136^*" 

+ ^ml -f Am2^ s^ge^'^'^-^cp) (nq) 
|mee| = \Jml + |Am§Jc^2C^3 

+ ^ml + Am2^ + I Ami J s\^cl^e^"' 
+m3S?3e"*(^-*cp) (10) (73) 

Obviously, these combinations strongly depend 
on the mass squared differences and mixing pa- 



rameters. In addition, there are two new indepen- 
dent phases a and /3, the Majorana phases. As a 
peculiarity, for the inverted ordering in Eq. (j73p . 
the first term dominates even for small values of 
TO3, which means that \mee\ cannot vanish. For 
the normal ordering in Eq. (l73l) . all the three 
terms could be of comparable magnitude. In this 
case, the three terms in Eq. (|73|) can be pic- 
tured as complex numbers to be added as vec- 
tors in the complex plane, where the angles be- 
tween the vectors are given by 2a and 2(/3 — i5cp) 
(see, e.g., Ref. [159] ). If the vectors are of similar 
length, the three vectors may add to zero, lead- 
ing to vanishing Oi>f5/3 decay. Therefore, an ob- 
servation of 0i>(3/3 decay is not guaranteed, even 
if the neutrinos are Majorana particles. We il- 
lustrate the functional dependence of \mee\ in 
Eq. ([73]) as a function of the lightest neutrino 
mass TOi (NO) or (10) in Fig. [131 In this 
figure, the uncertainties of phases and mixing an- 
gles within their current 3cr ranges are taken into 
account, which means that any point in the Oi//3/3 
regions could be allowed. For the NO, indeed 
0iyf3/3 may be suppressed, whereas for the 10, 
a non- vanishing \mee\ is guaranteed. For large 
TO, both cases are dominated by mi ~ TO3, i.e., 
|mee| increases linearly. Figures of this type are 
representative for presentations of 0vf3/3 decay; 
however, the details depend on the parameter as- 
sumptions. In the future, especially the combi- 
nation with the mass ordering measurements in 
long baseline experiments will be interesting. For 
example, the discovery of an inverted mass hierar- 
chy at long baseline experiments in combination 
with improved Oi^/3/3 bounds \mee\ ^ 0.01 eV ex- 
cluding the 10 solution would point towards Dirac 
neutrino masses. The current best bound |mee| ^ 
0.35 eV (90% CL) comes from the Heidelberg- 
Moscow collaboration [160j , where there is an un- 
certainty coming from the calculations of the nu- 
clear matrix elements (which we did not discuss 
here). An improvement of a factor of a few is 
expected, for example, from the GERDA exper- 
iment [163j . The required sensitivity to exclude 
the 10 solution might ultimatively be reached by 
the Majorana experiment |164j . 

Finally, indirect constraints to neutrino mass 
are nowadays obtained from cosmology. For ex- 
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ample, too much mass in relativistic neutrinos 
prohibits the growth of smaU scale matter per- 
turbations in the structure formation of the early 
universe, since energy is then transferred over 
longer scales. Cosmological tests of neutrino mass 
are typically sensitive to the sum of the neutrino 
masses m = J2i "^j ■ example, the gravita- 
tional effect of the neutrino masses does not de- 
pend on their mixings. The current limit depends 
on the data sets included, and is of the order 1 eV; 
see Ref. |165| for a very recent discussion. At 
this limit, neutrinos are quasi-degenerate, which 
means that mi ~ m2 ~ 7713 ~ 1/3 eV, which is 
used in Fig. 

7. Summary and conclusions 

In these lectures, we have discussed the phe- 
nomenology of the most important neutrino prop- 
erties, neutrino masses and mixings. In addition, 
we have shown a number of examples how to in- 
troduce new physics properties and how to test 
these in experiments. 

Most of the "canonical" neutrino properties 
are tested in neutrino oscillation experiments, 
namely the mixing angles, the Dirac CP phase, 
and the mass squared differences. We have shown 
that our current knowledge is based on two-flavor 
sub-sector measurements, often called the atmo- 
spheric oscillations, solar oscillations, and (short 
baseline) reactor oscillations. Of course, modern 
fits to data include the three-fiavorness of neutri- 
nos, but genuine three-flavor effects, such as lep- 
tonic CP violation, will be observable in the next 
generations of experiments. The key prerequisite 
for these measurements is the observation of ^13, 
see Ref. |166) for a recent review. Therefore, es- 
tablishing 9i3 > is the current top priority in 
neutrino oscillation physics. Beyond these canon- 
ical properties, future neutrino oscillation exper- 
iments could also probe new physics properties, 
such as coming from extra light sterile neutrinos 
or from LHC observable mediators. 

Apart from the mass squared differences, the 
absolute scale of the neutrino masses is needed 
to fix the mass spectrum. We have shown that 
different classes of experiments can probe the ab- 
solute neutrino mass scale, such as tritium end- 



point experiments, Qvj3p decay experiments, and 
cosmological tests. Whereas tritium endpoint ex- 
periments may be the most direct access to neu- 
trino mass, cosmological tests, which we have 
only shortly touched, provide strong indirect con- 
straints. Both of these tests work irrespective of 
the nature of neutrino mass, whereas the observa- 
tion of 0i>l3/3 decay will point towards Majorana 
mass terms. Majorana mass terms may be espe- 
cially interesting since they can be interpreted as 
lowest order perturbation of new physics. The 
simplest possibilities to obtain such Majorana 
mass terms are known as type I, II, or III seesaw 
mechanisms, depending on the heavy mediator 
leading to the low energy mass term. 

Except from neutrino sources such as the 
Earth's atmosphere, the Sun, or Earth-based ex- 
periments, we have also discussed cosmic accel- 
erators, such as extragalactic AGNs and GRBs. 
Although the main interest in these sources may 
be of astrophysical origin, they can be used to 
probe particle physics properties of the neutrinos 
as well. In particular, the observation of fiavor 
ratios in neutrino telescopes would help the ex- 
traction of such information. We have pointed 
out that in some extreme cases even information 
on neutrino flavor mixing may be obtained. It is, 
however, more likely that new physics properties, 
such as neutrino lifetime, will be probed using 
these sources. In order to illustrate the current 
state of research with respect to these measure- 
ments, he have discussed the sources, the propa- 
gation, and the detection, and we have sketched 
where current research requires further clarifica- 
tions. For example, the properties of the source 
change the flavor ratios as a function of energy, 
which, on the one hand, obscures the cleanliness 
of the flavor ratio, but, on the other hand, implies 
additional information. 

From the conceptual point of view, we have suc- 
cessively introduced neutrino oscillations in vac- 
uum, neutrino oscillations in constant matter, 
and then neutrino oscillations in varying mat- 
ter at the example of the Sun. Although we 
have only used the two-flavor limit, we have il- 
lustrated how the energy dependence of the solar 
oscillation effect can be reproduced in the sim- 
plest possible way. We have not included super- 
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nova neutrinos in the discussion, because collec- 
tive effects make a simple analytical treatment 
difficult. For the neutrino properties, we have in- 
cluded masses and mixings and some new physics 
examples. For other neutrino properties, such as 
electromagnetic dipole moments, see also lectures 
by Z.-z. Xing in this series p|. 

We conclude that the test of neutrino proper- 
ties has two interesting components. First of all, 
the leptonic mixing angles, the Dirac CP phase 
and the mass squared differences will be tested by 
future experiments, which are neutrino properties 
expected in all currently accepted scenarios. In 
addition, the absolute neutrino mass scale needs 
to be determined or constrained, maybe even so 
strong that the inverted ordering can be excluded. 
The process of 0iy(3(3 decay may also reveal the na- 
ture of neutrino mass, or constrain the contribu- 
tion of Majorana masses for the inverted order- 
ing. Apart from these guaranteed observations, 
neutrino physics may be one of the first places 
to search for new physics since massive neutrinos 
were not expected in the Standard Model. Es- 
pecially future precision oscillation experiments, 
such as neutrino factories, and the observation 
of high energetic neutrino fluxes at neutrino tele- 
scopes may help to constrain such new physics 
properties, either because of precision (neutrino 
factories) or because of the high energies (neu- 
trino telescopes). We therefore anticipate that 
neutrino physics could be good for further sur- 
prises. 
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